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y *************< 

/*  adtconv.c 

/*  adtconv.c  collects  data  from  GC  and  converts  it  to  digital  voltage 
/*  modification:  store  results  in  int  array  during  conversion  instead 
/*  of  storing  individual  byte  values 
/*  levels.  Base  address  is  set  for  0x20E. 

/*  created  4/8/85  B.  Lentz 

/* 

J*+1tir+*1Ht+1iiHHH**+++**+******+*****+********** 


7 

7 

7 

7 

V 

*/ 

V 


♦include  "bitset.h" 
♦include  <stdio.h> 

adtconv(dataval ,  gain, 
unsigned  int  *datava1; 
unsigned  int  gain; 
long  nconv; 
long  period; 

{ 

unsigned  char  temp; 
char  Hperiod; 
char  iperiod; 
char  Hno; 
char  Lno; 
long  loop: 
char  errorl; 
char  error?; 
unsigned  char  ad 1; 
unsigned  int  sum; 
char  setgain; 
float  volts; 
long  index; 
int  k; 

extern  unsigned  char 
extern  unsigned  char 


/*  bit  setting  macros  7 


nconv,  period) 

/*  raw  data  converted  to  decimal  7 
/*  value  used  in  voltage  calculation  7 
/*  number  of  data  conversions  to  perform  V 
/*  number  of  usee  between  points  (freq)  7 

/*  used  to  clear  out  register  7 
/*  high  byte  of  clock  period  7 
/*  low  period  of  clock  period  V 
/*  high  byte  of  #  of  conversions  7 
/*  low  byte  of  #  of  conversions  V 
/*  counter  for  conversion  loop  7 
/*  first  byte  of  error  code  7 
/*  second  byte  of  error  code  */ 

/*  low  byte  of  raw  data  */ 

/*  value  of  1st  pt  to  be  averaged  7 
/*  set  gain  on  board  7 
/*  voltage  of  data  value  7 


inportb(); 

outportb(); 


unsigned  char  start_chan  •  ' \0 1 ; 
unsigned  char  end_chan  *  ' \0 ‘ ; 


/*  only  using  one  channel  7 


if  (gain  »«  l) 
setgain  *  ' \0 ' ; 
else  if  (gain  **  10) 
setgain  *  '\1 * ; 
else  if  (gain  **  100) 
setgain  *  V ; 
else  if  (gain  **  500) 
setgain  *  ' \3 ’ ; 
else 

fprintf (stderr,  '\nILLEGAL  VALUE  FOR  GAIN\n  ); 
exit(0) ; 


if  ((STAT  REG  S  0x70)  !' 
I — { 


0) 


V.  O 
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fprintf (stderr,  "\nFATAL  ERROR- 1 1  legal  status  register  value\n"); 
fprintf (stderr,  "\nStatus  Register  value  is  %o\n“ ,  STAT_REG) ; 
exit(O) ; 


/*  stop  and  clear  */ 

C0MM_REG(CSTOP) ; 
temp  -  DATA  OUT; 
while(! ($TAT_REG  S  C0MM_WAIT)); 

COMMREG(CCLEAR)  ; 

whi  le(!  (STAT_REG  4  C0MM_WAIT)) ; 

COMM_REG(CCLOCK) ; 

Hperiod  *  (period/256)  4  0377; 

Lperiod  »  (period  -  (Hperiod  «  8) )  &  0377; 
whi le($TAT_REG  4  WRITE_WA1T) ; 

DATA  IN(Lperiod); 

whi  le(STAT_REG  4  WRITE  WAIT) : 

DATA_lN(Hperiod) ; 

whi  le(!(STAT_REG  4  COMM  WAIT))  ; 

COMM  REG(CSAD) ; 

whi  le(5TAT_REG  4  WRITE  WAIT) ; 

OATA  IN(setgain); 

whi  le(STAT_REG  4  WRITE  WAIT) ; 

DATAIN (start  chan); 

whi le(STAT_REG  4  WRITE_WA1T); 

OATA_IN(end_chan) ; 

Hno  *  (nconv/256)  4  0377; 

Lno  *  (nconv  -  (Hno  «  8))  &  0377; 
while  (ST  AT_REG  4  WRITEWAIT)  ; 

DATA  IN(Lno); 

while(STAT_REG  4  WRITE  WAIT) ; 

OATAJN(Hno)  ; 

while  (!  (STAT  REG  4  C0MM_WAIT)); 
COMM_REG(CRAD) ; 

index  -  nconv  »  3;  /*  divide  by  8  */ 

for (loop  -  0;  loop  <  index;  loop++) 

— { 

f or(k  »  0;  k  <  7;  k++) 
r-( 

I  while(! (STAT_REG  4  REA0_WA IT) ) ; 
adl  *  DATA_0UT; 

while(!(STAT_REG  4  PEA0_WAIT)); 
adl  «  OATAJXJT; 

whi le( ! (STAT_REG  4  REA0_WAIT) ) ; 
adl  *  DATA_0UT; 

whi  le( !  (STAT_REG  4  REAC_WAIT) ) ; 
dataval  [loop]"  *  (int) (0ATA_0UT  «  8)  +  adl; 
1—1 


'»  »**  j"*  /*  .'r  **  m  •  »  -  -  -  .  1  .  •>  *  . 
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whi le( ! (STAT_REG  &  C0MM_WAIT)); 
temp  *  $TAT_REG; 
if (temp  &  0x80) 
ioerr() ; 
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! ************************************************************************ i 


/*  analyz.c  */ 
/*  analyz.c  controls  sarple  collection  and  analysis  modules  */ 
/*  created  5/6/86  B.  Lentz  */ 
/*  V 
/*  modified  6/23/86  */ 


j *******************  ********************************************+********^ 

♦include  <$tdio.h> 

♦include  "parm.h" 

♦include  "targ.h" 

struct  sample  samp[MAXUNKS] ; 

struct  value  usr,  s3;  /*  contains  operating  parameters  */ 

main(angc,  angv) 
int  argc; 
char  *argv[); 

f 

int  frstrun;  /*  flag  first  run  of  day  */ 

printf ( "\nlnsert  data  diskette  into  Drive  B\n”); 
g_parm(&usr,  &s3); 
g_sspec(samp); 
frstrun  •  caldate(); 

if(atoi(argv[3])) 
usr.calflag  ■  TRUE; 
i  f  (usr.calf  lag) 

— { 

printf ('\nCALIBRAT ION  IS  REQUIRED\n\n") ; 

stdseq(); 

collect(l) ; 

— } 
else 

|  seqenc(frstrun)  ; 
j  collect(O); 

— } 


5 


/ 

♦ 

♦ 

♦ 

f 


;•  *.  \  >  a 


§  s 


■M 


i 


09-11-87  08:11:32  anltest.c 
Fri  09-11-87  10:04:51 


Pg  5 
of  205 
1-41 


j  ★  *■**★*•**★★  +  *■****★**★★**•*#******★**★★*★*******★■*★**■**•********#*********■**  j 

/*  analyz.c  controls  sample  collection  and  analysis  modules  */ 

/*  created  5/6/86  8.  Lentz  */ 

/*  modified  7/15/86  8.  Lentz  */ 

/*  printed  */ 

/****************************+*****************+************************* j 

#include  <stdio.h> 

♦include  "parm.h" 

♦include  "targ.h1' 


struct  sample  samp  [MAXIMS]  ; 

struct  value  usr,  s3;  /*  contains  operating  parameters  */ 


main(argc,  argv) 
int  argc; 
char  *argv[] ; 

-{ 

int  frstrun; 


/*  flag  first  run  of  day  */ 


printf ("\nlnsert  data  diskette  into  Drive  B\n"); 
g_parm(&usr,  Js3) ; 
gsspec(samp) ; 

/*  set  date  and  time  and  check  date  of  last  calibration  */ 
frstrun  «  caldateQ; 


if(atoi(argv[l])) 
usr.calflag  *  TRUE; 
if(usr.calflag) 


pri ntf ( "\nCAL IBRAT ION  IS  REQUlRED\n\n") ; 

stdseq(); 

collect(l); 

} 


/*  calibration  not  required  */ 
seqenc (frstrun) ; 
collect(O); 

} 


aldate.c 


caldate 


,  4*  +  **»****^**  +  #**•******★**★★***★*★★*  a#********************^ 

*/ 

1  c*  last  calibration  *, 

/%  3.  Lentz  *t 

* 

i 

>  * »♦**★***■********★*****★★★***★*********★★+********•***★*★★* 

t 


!"  month  of  last  calibration  */ 

/*  day  of  last  calibration  */ 

/*  calibration  date  in  days  since  previous  month  */ 
/ '  flap  first  run  of  day  */ 


*2tft*c%2d",  Smonth,  &day); 
.i,ete[0])  &&  (day  »*  ndate[l])) 

’•tSE; 


Time  %s\n",  usr.date,  usr.time); 
■Jibration  date  was  %s\n",  usr. caldate) ; 
c.?te.  '•%2d%*c%2d“ ,  Smooth,  Sday) ; 

-  nonth) 


i  -  da/;  >  7) 

!  ; 


i?[ij  +  mcn?riay[ndate[0]] ; 

i  >  /  ) 
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************************************* ********************************** **/ 


/ 

/*  coltest.c 

!*  collect  raw  data  and  analyze  it.  Display  retention  times  and  areas 
/*  test  response  to  trigger 
/*  2/6/86  B.  Lentz 
/*  modified  4/21/86 
/*  used  for  testing  with  temp  files 
/ 


V 

V 

V 

*/ 

V 

V 


^include  <5tdio.h> 

^include  "parm.h" 

♦include  "bitset.h" 

^define  MINMATCH  0.7 

col lect(calrun) 

int  calrun;  /*  is  this  a  calibration  run  ?  */ 


extern  struct  value  usr; 
int  position; 
int  npeaks; 
int  j; 

int  xmax [NPEAKS] ; 
float  area [NPEAKS] ; 
float  rettm; 
unsigned  int  *dataval; 
long  tarea[MAXSTDS] ; 
int  trettm[HAXSTDS] ; 
float  hit; 
float  match(); 

FILE  ‘pmptr; 

FILE  ‘floppy; 
char  rspns; 


/*  sample  position  */ 


/*  areas  for  all  peaks  */ 

/*  retention  times  for  all  peaks  */ 

/*  selected  areas  that  correspond  to  stds  &  target  */ 
/*  selected  retention  times  from  sample  V 
/*  match  factor  (pattern  ident)  */ 

/*  find  std  and  target  peaks  V 
/*  used  to  send  data  to  printer  */ 

/*  used  to  send  data  to  floppy  file  */ 

/*  first  letter  of  user  response  */ 


/*  open  printer  file  */ 
if((pmptr  «  fopen( "PRN:",  "wb"))  ■*  0) 

■-{ 

fprintf(stderr,  "\nError  transmitting  data  to  printer  \n"); 
return; 

-} 

/*  send  carriage  return  and  form  feed  signal  to  Epson  printer  */ 
fprintffpmptr,  "%c%c",  12,  13); 

print f("\nLoad  Blank  in  position  1  and  Start  Tekmar  Purge  and  Trap\n'); 

l — { 

printf ("\n\nls  the  'Purge  Complete'  light  illuminated  on  the  LSC?  (y  or  n)  :"); 
[  while(! (rspns  -  replyO)); 
if(rspns  <  0)  /*  exit  V 

return; 

usr. invalid  *  FALSE; 

Composition  *  0;  position  <  usr.sampcnt;  position*-*) 

-{ 

setdate() ; 

printf ("\n\n\tAnalyzing  SAMPLE  f  %d\t%s\n",  position-*-! ,  usr.mwber[position]) ; 
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ptisuion,  >uhx,  area); 
position) ; 


"iiv-  \n\tpeak\tretention  time\t\tarea\n\n") ; 
npeaks;  j++) 

■  f;  1  /  !'isr.freq*60.0) ; 

;"rf  \t  %d  \t  %.3f  \t\t%.0f\n",  j,  rettm,  area[j]); 


npeaks,  xmax,  area,  0); 
it  ion,  1 ,  -1.0,  0.0) ; 

rr,  r  k‘'*  WARNING  ***  chromatography  problem  \n  Detected  only  %d  peaks\n',  npeaks); 
,n'*‘  WARNING  ***  chromatography  problem  \n  Detected  only  %d  peaks\n",  npeaks); 

• ;  -  vj) 

•osu’-jO.  rpeaks,  xmax,  area,  0); 

position,  1,  -1.0,  0.0); 

•  ro-r,  \r,**»  WARNING  ***  Sample  is  too  complex  for  automated  analysis\n"); 

-r-r  \n**»  WARNING  ***  Sample  is  too  complex  for  automated  analysis\n”); 


•jilnpcaks,  xmax,  area,  tarea,  trettm); 

*  >cn.  npeaks,  xmax,  area,  hit); 
r.  "  nstMatch  factor  %f\n" ,  hit); 

‘•HATCH) 

■•r ir ,  \n\n***  WARNING  ***  Problem  with  chromatography  -  Internal  Standard  peaks  not  found\n"); 
c.v-rr,  ‘\n\n***  WARNING  ***  Problem  with  chromatography  -  Internal  Standard  peaks  not  found\n"); 

■ 1  -  twi;e; 

;,r r.  position,  trettm,  tarea,  icalrun); 

■  cea*s  V 

::■■■  ‘  it'Cn  */ 

'  turn  and  form  feed  signal  to  Epson  printer  */ 

12,  13); 
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/ft***********************************************************************  j 

/*  eollect.c  */ 
/*  collect  raw  data  and  analyze  it.  Oisplay  retention  times  and  areas  V 
/*  test  response  to  trigger  */ 
/*  created  2/6/86  B.  Lentz  */ 


7*  modified  4/21/86  V 

/+++*+*#**#******************+*********************♦#********************/ 


include  <stdio.h> 

♦include  "parm.h" 

♦include  "bitset.h" 

collect(calrvn) 
int  cal  run; 

-{ 

extern  struct  value  usr; 
int  position;  /*  sa 

int  npeaks; 
int  j; 

int  xmax[NPEAKS]; 
float  area[NPEAKS] ; 
float  rettm; 
unsigned  int  *dataval; 
long  tarea[MAXSTDS] ;  /*  s t 

int  trettm[MAXSTDS] ;  /*  s( 

float  hit;  /*  ma 

float  match();  /*  fi 

FILE  *pmptr;  /*  us 

FILE  *floppy;  /*  us 

char  rspns;  /*  fi 

unsigned  int  core;  /*  an 

extern  unsigned  int  coreleft(); 


sample  position  */ 


/*  selected  areas  that  correspond  to  stds  &  target  */ 
/*  selected  retention  times  from  sample  */ 

/*  match  factor  (pattern  ident)  */ 

/*  find  std  and  target  peaks  */ 

/*  used  to  send  data  to  printer  */ 

/*  used  to  send  data  to  floppy  file  */ 

/*  first  letter  of  user  response  */ 

/*  amount  of  memory  available  */ 


if((pmptr  -  fopen("PRN:",  "wb"))  *■  0) 

p{ 

fprintf (stderr,  "\nError  transmitting  data  to  printer  \n" ) ; 
return; 

fprintf (pmptr,  "%c%c\  12,  13); 

printf ("\nLoad  Blank  in  position  1  and  Start  Tekmar  Purge  and  Trap\n  ); 


— do  { 

printf ("\n\nls  the  'Purge  Complete'  light 
— }  while(! (rspns  *  replyQ)); 


illuminated  on  the  LSC?  (y  or  n) 


if  (rspns  <  0) 
return ; 


initdio() ; 
iffusr. verbose) 

-{ 

printf("\ngain  %u  nconv  %0  period  %0\n",  usr. gain,  usr.nconv,  usr.period); 
printf ("runtime  Vf  min\tfrequency  Vf  pts/sec\n\n" ,  usr. runtime,  usr.freq); 

printf ("threshold  values  Vd\t  %d\t  Vd\t  %d\n",  usr.thresh(l] ,  usr. thresh [2] ,  usr.thresh(3],  usr.thresh{4]) ; 
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core  *  coreleft() ; 
if(core  <  (usr.numpt  *  2)) 

— { 

fprintf(stderr,  "\n\nFATAL  ERROR  not  enough  memory  available  for  data  arrays\n") 
return; 

usr. invalid  -  FALSE; 

for(po$ition  *  0;  position  <  usr.sanpcnt;  position*-*-) 


setdate() ; 

printf( "\n\n\tAnalyzing  Sample  #  %f\t%s\n“ ,  position+1,  usr.mmber[position]) ; 
if  ((dataval  =  alloc(usr.nunpt  *  2))  *«  0) 


fprintf (stderr,  "not  enough  memory  available\n") ; 

return; 


dataval  *  Sdatava![0]; 


setdio(01,  CDKXXfT) ; 

setdio(0,  XTRIG);  /*  wait  for  trigger  */ 
printf ("\nCollecting  Data  \n"); 

adtconv (dataval,  usr. gain,  usr.nconv,  usr. period); 

smooth  (dataval ,  usr.nunpt,  usr. smooth); 

npeaks  *  pk  det(usr.mjnpt,  dataval,  Susr,  xmax,  area); 
pmres(pmptr,  position); 
if(npeaks  >  0) 


fprintf (pmptr,  "\n\tpeak\tretention  time\t\tarea\n\n"); 
for(j  *  1;  j  o  npeaks;  j-*-*-) 

-{ 

rettro  *  xmax[j]/(usr.freq*60.0) ; 

fprintf(pmptr,  "\t  %d  \t  %.3f  \t\t%.0f\n”,  j,  rettm,  area[j]); 


if(npeaks  <  3) 

r< 

stof lpy(position,  npeaks,  xmax,  area,  0); 

fprintf  (pmptr,  "\n***  WARNING  ***  chromatography  problem  \n  Detected  only  %d  peaks\n",  npeaks); 

fprintf (stderr,  "\n***  WARNING  ***  chromatography  problem  \n  Detected  only  %d  peaks\n’,  npeaks); 

else  if (npeaks  >  (NPEAKS-1)) 

r"f 

stoflpv(position,  npeaks,  xmax,  area,  0); 

j  fprintf f pmptr,  "\n**‘  WARNING  ***  Sample  is  too  complex  for  automated  analysis  n  ); 

j  ;  fprin'f. 'stderr,  \n***  WARNING  ***  Sample  is  too  complex  for  automated  analysis  n  ); 


|  hit  '  raatch(npeaks,  xmax,  area,  tarea,  trettm); 
I  stot ipyfposition,  npeaks,  xmax,  area,  hit); 
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fprintf(pmptr,  "\n\tMatch  factor  %f\n",  hit); 
if (hit  <  0.7) 

r{ 

fpnntf(pmptr,  "\n\n*»*  WARNING  ***  Problem  with  chromatography  -  Internal  Standard  peaks  not  fordin'1 
fpnntf(stderr,  "\n\n***  WARNING  ***  Problem  with  chromatography  -  Internal  Standard  peaks  not  ford\n" 
usr. invalid  •  TRUE; 

else 

sampler(pmptr,  position,  trettm,  tarea,  Jcalrun); 
free(dataval); 

fprintf(pmptr,  "%c%c",  12,  13); 
fclose(pmptr); 

-} 
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*+-»******+*****★******-***  ***************************************  ******** 


/ 


V 
*/ 
*/ 

V 

revised  6/9/86  */ 

****** **************************************************************** j 


/*  chgtgt. c 

/*  chgtgt.c  Change  target  paramters 
/*  created  8/26/85  8.  Lentz 


♦include  <stdio.h> 
♦  •include  'parm.h" 
♦include  targ.h 


crgtgtO 


extern  struct  value  usr;  /*  user  parameters  */ 
struct  saddle  samp[MAXUNKS] ;  /*  target  paramters  */ 

/*  selected  option  to  change  */ 

/*  nunerical  value  of  character  input  */ 
/*  selected  target  cmpd  to  change  */ 


cha 

lorn 


cpi.'on; 

Inum; 


targ; 

"y  a*-oi(), 


;  ong 

float  fnutrfcer; 
float  q  tnum(); 
chai  string[30]; 
*nt  i; 

char  newval  [20] ; 
char  .name [20]; 
FILL  *f i  1  eptr; 


/ 


/ 


*  floating  point  number  */ 

/*  get  floating  point  number  V 
user  input  */ 

/*  loop  counter  */ 

/*  new  value  of  parameter  */ 

/*  parameter  name  */ 


/*  clear  the  screen  */ 

system("cls"); 

printf]  Reading  current  values  of  parameters\nPlease  wait...\n"' 

yjspec(sanp) ; 
tarq  -  1; 


f'tfileptr  *  fopen(targfile[targ] ,  ”rw” ) )  *■  0) 


fprintffstderr,  "error  opening  parameter  fi le\n") ; 

return; 


whi  lu(0  *«  0) 


I  ! 


system] 1 
pnntf  (" 
printf (" 
printf  (' 
printf (" 
print*] ' 
printf (" 
pr intf ( " 
priltf] 
printf (" 
iinnif  {" 
printf  (' 
pr'ntf  ]' 


els"); 

\n\n\na)  Target  compound  is  %s\n" 


samp [targ] .name) ; 


60.0)); 

samp[targ] .cutoff) ; 


Retention  Time  is  %.2f  min\n",  ]samp[targ] ,rt)/(usr.freq 
Units  of  Concentration  %s\n",  samp[targ] .units); 

Haxiirun  allowable  water  sample  limit(waminq  limit)  is  %.2f\n" 

Spike  concentration  is  %.2f  \n",  sampftargj. spike); 

Concentration  of  ’ss  in  standard  samples:\n ",  $amp[tang].name); 

\tConcentration  in  standard  A  is  %.2f\n",  samp[targ] ,conc[l]) ; 

\tConcentration  in  standard  B  is  %.2f\n",  sarp[targ] .conc[2]) ; 

\tConcentration  in  standard  C  is  %.2f\n",  samp[targ].conc[3]); 

g)  Precision:  maximum  difference  between  duplicate  spike  samples  ■  %.2f\n",  samp [targ] .precisn) ; 

h)  Accuracy:  maximum  difference  between  actual  and  measured  results  •  %.2f\n",  samp  [targ]. aceurcy); 

i)  Top  concentration  of  linear  range  *  %.2f  %s\n",  sampftarg] .dilution,  samp[targ]  .uiits) ; 
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printf("j)  Quantisation  limit  is  %.lf  %s\n",  samp[targ].qtlimit,  sanp[targ]  .units) ; 
printf(“x)  Exi t\n" ) ; 

printf ("\n£nter  letter  of  option  to  change:  "); 
option  -  getchar(); 
get  char  ()  ; 

switch(option) 

r-1 

case  'a' : 
case  'A': 

r{ 

printf ("\nChanging  name  of  target  compound  #%0\n" ,  targ); 
printf ("\nEnter  new  compound  name:  "); 
if (fgets (string,  25,  stdin)  **  0) 

r{ 

printf ( "\nlnvalid  compound  name\nTarget  name  not  changed\n"); 
break; 

L-} 

wrparm("name ",  string,  fileptr); 
break ; 

M 

case  ‘b1 : 
case  'B' : 

r{ 

printf{"\n\n  Changing  Retention  Time  (in  minutes)  of  target  compound  %D",  targ); 
if((fnurber  *  gfnum(string))  <  0.0) 

r{ 

printf ("\nRetention  time  of  %s  not  changed\n",  samp(targ] .name) ; 
break; 

L-} 

inum  *  fnunber  *  usr.freq  *  60.0; 
if  (fnunber  >  usr. runtime) 

— { 

printf ("\nlnval id  Retention  Time  -  Larger  than  run  time\n'); 
printf ("\nRetention  time  of  %s  not  changed\n",  samp[targ].name) ; 
break; 

sprintf(newval,  "%d",  inum); 
wrparm("rettime" ,  newval,  fileptr); 
break; 

L) 

case  ‘ c ' : 
case  'C' : 

r{ 

printf ("\nChanging  Units  of  Concentration^ ')  ; 
printf ("\n£nter  new  unit  name:  "); 
if(fgets(strinq,  25,  stdin)  *■  0) 

printf ("\nlnval id  name,  unit  name  not  changed.n  ); 

break ; 

L-} 

wrparm( "name",  string,  fileptr); 

break; 

L} 

case  'd ' : 
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case  'O': 

r{ 

printf ("\n\n  Changing  Maximum  allowable  limit  for  concentration"); 
if((inuit  »  getnum(newval))  <  0) 

r{ 

printf("\n  Invalid  cutoff  limit,  value  not  changed\n"); 
break; 

wrparm("cutoff",  newval,  fileptr); 
break; 

L-} 

case  ' e* : 
case  'E': 

p{ 

printf("\n\n  Changing  concentration  of  spike  added1’); 
if((imjn  -  getm»(newval))  <  0) 

r{ 

printf("\n  Invalid  spike  concentration,  value  not  changed\n"); 
break; 

L-} 

wrparm("spike'‘,  newval,  fileptr); 
break; 

case  'f ' : 
case  1 F’ : 

r{ 

printf ("\n\n  Changing  concentration  of  %s  in  standard  A",  samp[targ].name) ; 
if((inun  -  getnun(newval))  <  0) 

r{ 

printf ("\n  Invalid  concentration,  value  not  chenged\n” ) ; 
break; 

L-} 

wrparm("stdconcl",  newval,  fileptr); 

printf ("\n\n  Changing  concentration  of  %s  in  standard  B",  samp[targ].name); 
if((inun  *  getnum(newval))  <  0) 

H 

printf ("\n  Invalid  concentration,  value  not  changed\n") ; 
break; 

wrparm("stdconc2" ,  newval,  fileptr); 

printf ("\n\n  Changing  concentration  of  %s  in  standard  C",  samoTtarq] .name) ; 
if((inu»  ■  getnum(newval))  <  0) 

r{ 

printf ("\n  Invalid  concentration,  value  not  changed\n"); 
break; 

wrparm("stdconc3" ,  newval,  fileptr); 
break ; 


case  'g‘: 
case  'G' : 
r-{ 

printf("\n\n  Changing  precision  control  value'"); 
if((inu«  •  getnum(newval))  <  0) 

! — { 
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printf("\n  Invalid  precision,  value  not  changed\n"); 
break; 

M 

wrparm("precisn",  newval,  fileptr); 
break ; 

*-} 

case  1 h 1 : 
case  'H': 

printf("\n\n  Changing  accuracy  control  value"); 
if((inun  »  getnum(newval))  <  0) 
r-{ 

printf("\n  Invalid  accuracy,  value  not  changed\n"); 
break; 

L-} 

wrparm("accurcy",  newval,  fileptr); 
break; 

M 

case  'i' : 
case  ' I' : 
r-{ 

printf("\n\n  Changing  top  of  linear  range.  Concentration  should  be  in  %s",  sampftarg). units); 
if((inun  -  getnum(newval) )  <  0) 

printf ("\n  Invalid  concentration,  linear  range  not  changed\n“); 
break; 

L} 

wrparm("dilution" ,  newval,  fileptr); 
break ; 

U 

case  ' j 1 : 
case  ' J' : 

printf ('\n\n  Changing  guantitation  limit.  Concentration  should  be  in  %s",  samp[targ], units); 
if ( (inun  -  getnum(newval))  <  0) 

— { 

printf ("\n  Invalid  concentration,  guantitation  limit  not  changed\n"); 
break; 

wrparm("qtlimit" ,  newval,  fileptr); 
break ; 


case  ' x ' : 
case  ' X 1 : 
r-{ 

fclose(fileptr); 
return; 
break ; 

M 

default: 

p{ 

printf ("\n\tlnval  id  selection\n") ; 
break ; 

M 
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j **<•******»*♦************************************************* 

/*  change. c 

/*  change. c  Charges  values  of  user  parameters  in  "parm1  file 
/*  created  8/15/85  B.  Lentz 

/* 

/*  revised  6/12/85  ^ 

♦include  «3tdio.h> 

♦include  "parm.h" 

struct  value  usr,  s3; 


main() 

-{ 

char  newval [10] ; 
char  string  [10]; 
int  base  *  0; 
long  inumber; 
long  atoi () ; 
float  f number; 
float  g_fnum(); 
long  getram(); 
long  i ; 

char  name [20]; 
int  option; 
int  level; 
int  rspn; 

FILE  *fptr; 


/*  level  of  option  to  change  */ 


if((fptr  *  fopen(targfile[0] ,  "rw"))  **  0) 

fprintf (stderr,  "\tERR0R  opening  parameter  file\n"); 
return; 

-} 

while(0  **  0) 

r{ 

g_parm(&usr,  is3); 
pmparm(&usr)  ; 

printf("\nEnter  letter  of  option  to  charge:  "); 
option  ■  getcharQ ; 
getchar() ; 

switch(option) 

H 


case  'A' 


printf ("\h\n\n") ; 

printf("\n\nChanging  Verbose  0ption\n'); 
do 

printf{"\nSet  Verbose  option  ON  ?  (y  or  n):  '); 
rspn  •  getchar(); 
getchar() ; 


main 
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if(rspn  **  'y'  |  rspn  **  *  V ) 

s  trxpy(  newval ,  "ON"); 
else 

strcpy(newval ,  "OFF"); 

— }  while  (rspn  <  1); 
wrparm("verbose" ,  newval ,  fptr) ; 
break; 

case  'b': 
case  ' B‘ : 
r-{ 

printf ("\n\n\n") ; 

printf ("\n\nChanging  chromatographic  sample  time\n"); 
if((getnum(newval))  <  0) 

— { 

printf ("\nChromatographic  sample  time  not  changed\n"); 
break; 

wrparm(" runtime",  newval,  fptr); 
break; 

L} 

case  ' c ' : 
case  'C' : 
r~{ 

printf ("\n\n\n"); 
for (i  ■  1;  i  <  4;  i++) 

r{ 

sprintf(name,  "stdRT%d",  i); 

printf ("\nChanging  Retention  Time{in  minutes)  of  internal  standard  f%d",  i); 
if((fniwber  *  g_fnim(string))  <*  0.0) 

— { 

printf ("\nRetention  time  of  internal  standard  #%d  not  changed\n",  i); 
break; 

L_} 

inumber  •  fnunber  *  usr.freq  *  60.0; 
if  (fnumber  >  usr. runtime) 

r{ 

printf (" Invalid  Retention  Time  -  Larger  than  run  time\n"); 
break; 

L} 

sprintf (newval,  "%d",  inimber); 
wrpanm(name,  newval,  fptr); 

L-} 

break ; 

L-} 

case  ’ d 1 : 
case  'D': 

_{ 

printf ("\n\n\n") ; 
for ( i  •  1;  i  <  4;  i++) 

printf( "\nChanging  AREA  of  internal  standard  »%d  (for  match  calculation)  ,  i); 
sprintf (name,  "stdarea%d",  i); 
if((getnum(newval))  <*  0) 
i — { 

19 
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printf("\nArea  of  standard  #%d  not  changed\n",  i); 
break ; 

L-} 

wrparm(name,  newval,  fptr); 

break; 

*-} 

case  ' e '  : 
case  'E' : 

printf( '\n\n\n") ; 

printf (”\n\nChanging  minimun  area  required  for  standard\n") ; 
if ((getnum(newval))  <■  0) 
r-{ 

printf ("\nMininun  area  required  of  standards  not  changed\n"); 
break ; 

>-} 

wrparm( "minarea" ,  newval,  fptr); 
break; 

case  ' t ' : 
case  'T': 

— { 

chgtgt(); 
break ; 

case  ’ x ' : 
case  'X' : 
r-{ 

fclose(fptr) ; 
return; 
break ; 

default: 


1  printf( 


\n\t Inval id  selection\n ') ; 
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j ******************************* ******************** ************** ******* j 


/*  colunn.c  */ 
/*  Calculates  and  prints  colunr  factor.  Monitoring  this  factor  gives  7 
/*  an  indication  of  the  condition  of  the  colunr.  The  value  decreases  */ 
/*  as  the  colunr  deteriorates.  */ 
/*  created  8/28/86  B.  Lentz  */ 
/*  V 
/*  Calling  Modules:  collect. c  V 
/*  Modules  Called:  none  V 
/*  V 


/******★**★****★********★*★★*★*★*★*★★*★***★★**★★★**★****★★*★*★****★*****★ 


♦include  <5tdio.h> 


coluir(pmptr,  ymax,  tarea,  stindx,  nstd) 


FILE  *pmptr; 
int  ymax[]; 
long  tarea[]  ; 
int  stindxQ ; 
int  nstd; 


/*  used  to  send  data  to  printer  */ 

/*  peak  heights  */ 

/*  selected  areas  that  correspond  to  stds  &  target  V 
/*  array  index  to  standard  peaks  */ 

/*  number  of  internal  standards  */ 


r~{ 
int  i ; 

float  cfactor;  /*  column  factor  */ 


fprintf(pmptr,  "\n\tColunr  factors:"); 
for(i  •  1;  i  <=■  nstd;  i++) 

r{ 

cfactor  «  (float)  ymax[stindx[i]]/tarea[i] ; 
fprintf(pmptr,  "\t%.4f",  cfactor); 

M 

■— } 
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ft*********** ******** *★**★**★***★**********★**★************************/ 


/*  copsys.c 

/*  format  a  new  diskette  and  copy  operating  system  and 
/*  datasystem  files  to  it 
/*  created  7/21/86  B. Lentz 


/* 


/ 


***+★★*** 


♦include  <stdio.h> 
♦define  SPACE  '  ' 


main() 

-{ 

int  ntime[4]; 
int  ndate[4] ; 
char  string [100]; 


printf ("\n\n\nThis  utility  creates  a  new  system  operating  disk\n"); 
printf ("These  disks  should  be  replaced  month ly\n“); 
printf ("\n\nlnsert  a  NEW  diskette  into  Drive  B\nReady  (y  or  n):  "); 
while(!reply()) 

printf ("New  Diskette  in  Drive  B?  (y  or  n):  ”); 
printf (“\nPlease  Wait..."); 
system("type  response  |  format  b:/s  >NUL"); 
system("copy  a:*.*  b:  >  NUL"); 
system("del  b:*.dat  >  NUL"); 
sysdate(ndate,  ntime); 

sprintf (string,  "copya.QC%d%d.dat  b : " ,  ndate[0] ,  ndate[2}); 
system(string) ; 

printf("\n\nNew  system  diskette  is  ready  for  use  in  Drive  A\n"); 
printf ("\nRemove  diskette  from  Drive  B  and  place  a  label  on  it\n"); 
printf ("Place  QC  archive  diskette  in  Drive  B\n\nReady  (y  or  n):  "); 
whi le(! reply ()) 

printf ("QC  Archive  Diskette  in  Drive  B?  (y  or  n):  "); 
printf ("\nPlease  Wait..."); 
system("copy  a:*.dat  b:  "); 

printf ("\nRemove  QC  archive  diskette  from  Drive  B  and  place  data  diskette  in  Drive  B\n") ; 


printf ("\nReady  to  continue?  (y  or  n): 
while(!reply()) 

printf ("Ready  to  continue?  (y  or  n): 


'); 
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j *************************** ************************************** *^ 

/*  error.c  */ 

/*  print  error  messages  */ 

/*  created  5/9/86  B.  Lentz  V 

/***************♦**♦**********************♦************************/ 

#include  <stdio.h> 

error(code) 
int  code; 

switch (code) 
r-{ 

case  0: 

r{ 

fprintf(stderr,  "\n\nProblem  with  chromatography-  Standard  peaks  not  found\n"); 
break; 

*-} 
case  1: 
default: 

— ^ 

fprintf(stderr,  "\n\nOata  System  Error\n"); 
break; 

L> 

M 

L> 
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/*  filelist.c  */ 

/*  get  list  of  data  files  on  floppy  disk  */ 

/*  created  7/1/86  B. Lentz  */ 

7 

/**»**************»*****«************t,***i,-n,in,tittin,**1,i,t»tin,tin,**ini,in,**tj 


#include  <stdio.h> 

♦define  TRUE  (0—0) 

♦define  FALSE  ! TRUE 

filelist() 

-{ 

extern  char  *filedir();  /*  returns  list  of  matching  filenames  */ 
char  ‘list;  /*  list  of  filenames  */ 

char  ‘next; 

char  f i 1 espec [20] ;  /*  file  specification  */ 

char  filename [20]; 

int  cot nt;  /*  filename  count  used  to  divide  output  in  thirds  */ 

int  length; 

FILE  ‘pmptr; 

int  reply();  /*  get  yes  or  no  reply  from  user  */ 


/*  get  yes  or  no  reply  from  user  */ 


printf(''\n\nlnsert  Data  diskette  into  Drive  B\nReady?  (y  or  n):  “); 
while(!reply()) 

printf("Oata  Diskette  in  Orive  B?  (y  or  n);  "); 
printf ("\n\n\tGathering  filenames  from  disk  B\n\tP1ease  Wai t\n" ) ; 
system("ls  -m  b:*.*"); 

printf ("\nDo  you  wish  to  have  a  printout  of  the  available  files?  (y  or  n):  '); 
if (reply()) 

— { 

if((pmptr  »  fopen("PRN:",  "w"))  »»  0) 
r-{ 

fprintf(stderr,  "\nError  transmitting  data  to  printer  \n"); 
return; 

strcpy(filespec,  "b:*.*"); 

if ( (1 ist  *  filedir(fi lespec,  0))  ■»  NULL) 

— ^ 

printf("\nThere  are  no  data  files  archived  on  the  floppy  disk\n") ; 
return; 

count  *  0; 

for(next  *  list;  ‘next  !*  NULL;) 

r{ 

fprintf (pmptr,  "%s\t\t\t",  next); 
next  +-  strlen(next)  +  1; 


if(++count  >  1) 

-{ 

fprintf  (pmptr,  "\n"); 
count  *  0; 

-i 


L> 

fc  lose  (pmptr)  ; 
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j _ free(list); 


printf("\nENTIR  filename  to  retrieve:  "); 
fgets (filename,  19,  stdin); 
while(stranp(fil«’»ie.  "\n")  !*  °) 


length  ■  strlen(filename); 
filename[length-l]  *  '\0‘ ; 
retriev (filename); 

printf ("\nENTER  filename  to  retrieve:  ”); 
fgets (filename,  19,  stdin); 

-} 


$  $ 
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^***+******+*************************************************************  j 

/*  fixseq.c  */ 

/*  correct  any  errors  made  during  sequence  entry(entering  sample  */ 

/*  numbers  and  positions  for  the  purge  and  trap)  */ 

/*  created  4/4/86  B.  Lentz  revised  6/11/86  */ 

/*  V 

j  to********************************** ************************************ / 


♦include  <stdio.h> 
♦include  "parm.h" 


fixseq(frstpos,  lastpos) 

int  frstpos;  /*  first  position  that  can  be  changed  */ 

int  lastpos;  /*  last  sample  position  on  autosampler  to  be  used  V 

-{ 

extern  struct  value  usr; 

extern  char  label Q;  /*  date  and  time  information  used  in  name  */ 


int  pos;  /*  sample  position  on  purge  and  trap  */ 

int  type;  /*  sample  type  */ 

int  num;  /*  temporary  storage  for  int  */ 

char  rspns;  /*  first  letter  of  user  response  to  query  7 

char  string  [20];  /*  temporal  string  storage  */ 

char  filename [10] ; 


system("cls"); 

strcpy (filename,  "CON:"); 

pmseq(filename); 

printf ("\n\n\t  CHANGE  ANY  SAMPLE  NUMBERS?  (y  or  n):  "); 
if(rspns  *  replyO) 

— { 

p-do  { 

printf ("\n\tEnter  POSITION  NUMBER  of  sample  to  be  changed:  "); 
pos  »  select(); 


— }  while(pos  <  frstpos  ||  pos  >  lastpos); 
pos-; 

if (usr.type[pos]  *•  SAMPLE) 

-{ 

sprintf(string,  "del  b:%s",  usr.number[pos]); 
system(string) ; 

-} 

printf(”\n\t\t\tSample  Type  Numbers\n\t  0)  BLANK\t\t\t4)  WATER  SAMPLE\n"); 

printf("\t  1)  STAN0AR0  #A\t\t5)  SPIKE\n") ; 

printf ("\t  2)  STANDARD  #B\n\t  3)  STANDARD  ♦C\n\n"); 

printf("\tPosition  %d\tENTER  sample  TYPE  number:  ",  pos+1); 

type  *  select{); 

while(type  <  0  J type  >  5) 

-{ 

printf ("\t\tENTER  sanple  TYPE  number:  ",  pos+1); 
type  *  select(); 

' — } 

samptype(Spos,  type); 

-} 


■^<0-  - 


-j-- ^ 
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2  /*  fmtdis.c  */ 

3  !*  format  a  new  diskette  for  data  collection  */ 

4  /*  created  8/8/86  8. Lentz  V 

5  /*  */ 


^  j ******■**♦★*********★%★*******★*******★****★******★*****★****♦* ********** j 

7 

8  linclude  <stdio.h> 

9 

10  main() 

11  H 

12  printf ( "\n\n***WARNIHG***  All  data  on  diskette  in  drive  B  will  be  lost\n“); 

13  printf ("\n\nlnsert  a  new  diskette  into  Drive  B\nReady?  (y  or  n):  ”); 

14  whi1e(!reply()) 

15  printf ("New  Diskette  in  Drive  B?  (y  or  n):  "); 

16  printf( "\nPlease  Wait..."); 

17  system("type  response  |  format  b:  >  NUL"' ) ; 

18  printf("\n\nNew  data  diskette  is  ready  for  use  in  Drive  B\n"); 

19  print f("\nReady  to  continue?  (y  or  n):  "); 

20  wtiile(!  reply  ()) 

21  printf ("Ready  to  continue?  (y  or  n):  "); 

22  1— } 

23 
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y  imHt*********************************************************************} 

2  /‘get  current  system  date  and  time 

3  /*  created  8/7/85  revised  *' 

4  /*  B.  Lentz  *' 


6 

7  getdate(date,  time) 

8  int  dateQ; 

9  int  time[]; 

10  r{ 

11  struct  regval  {  int  ax,bx,cx,dx,si,di,ds,es;  }  srv; 

12  mtndat[4]; 

13 

14  srv. ax  *  0x2a00; 

15  sysint21(&srv,&srv) ; 

16  date [2]  •  srv.cx;  /*  year  V 

17  ndat [1]  ■  srv.dx; 

18 

19  srv. ax  -  0x2c00; 

20  svsint21(&srv,&srv); 

21  ndat[2]  *  srv.cx; 

22 

23  date[0]  •  ndat [1]  »  8;  /*  month  */ 

24  date[l]  •  ndat[l]  fi  Oxff;  /*  day  */ 

25  time[0]  -  ndat[2]  »  8;  /*  hour  */ 

26  time[l]  •  ndat [2]  &  Oxff;  /*  minutes  */ 


->  J 
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I  mt  *******  *************************************************/ 

/*  getname.c  */ 

/*  check  for  valid  filename  -  alphanumerics  only  */ 

/*  created  8/13/86  B.  Lentz  */ 

j +++*++*+****+ it***************************** ******* *++***++*/ 

♦include  <stdio.h> 

getname(filename) 
char  ‘filename; 

-{ 

char  ‘string,  *ptr; 
char  c; 

FILE  *fptr; 

extern  char  *calloc(); 


string  *  calloc(20,  1); 
ptr  *  string; 
if(! (gets(filename,  16))) 
pri ntf ( "\n\t\tERR0R  reading  sample  number:\t  %s 
strcpy(string,  filename); 


filename) ; 


while((c  *  ‘string)  !»  NULL) 

-{ 

if(isalnum(c)  **  0) 

printf ("\n\t\t\tlnval id  character  in  %s\n",  filename); 
printf("\t\t\tRE-ENTER  sample  IDENTIFICATION:  "); 
free(ptr) ; 
retum(l); 

*— •} 

string -w-; 

I-} 

string  *  ptr; 

strcat(filename,  ".dat"); 
sprintf(string,  "b:%s",  filename); 
if((fptr  •  fopen(string,  r"))  !*  0) 

r{ 

fclose(fptr) ; 

printf ("\n\t This  filename  already  exists  on  disk,  please  select  another  name'n"); 
printf ("\t\t\tRE-ENTER  sample  IDENTIFICATION:  ); 
free(ptr) ; 
retum(l) ; 

L} 

else 

r-{ 

if((fptr  »  fopen(string,  "w"))  **  0) 

— { 

/*  out  of  space  on  floppy  V 
free(ptr); 
retum(2) ; 

fclose(fptr); 
free(ptr) ; 
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J  j  ***********#******************+*****************  it**********************/ 

2  /*  getnun.c  */ 

3  !*  getnum.c  Pronpt  user  for  nunerical  input  and  then  verify  it.  */ 

4  /*  If  valid  string  is  found,  returns  length  of  number  string,  */ 

5  /*  otherwise  returns  zero.  */ 

6  /*  created  8/15/85  B.  Lentz  */ 

7  /*  V 

g  j a*************************************************************** ******** j 

9  #include  «stdio.h> 

10 

11  long  getnun(newval) 

12  char  ‘newval; 

13  r{ 

14  int  i;  /*  counter  for  loop  */ 

15  long  nun;  /*  numerical  value  of  parameter  V 

16  long  atoi () ; 

17  char  string[20];  /*  input  string  containing  new  value  */ 


printf  ("\n\tEnter  new  value:  '') 
if(fgets(string,  20,  stdin)  **  0 
printf ("\nerror  in  input\nM); 
sscanf(string,  "its" ,  newval); 
nun  «  atoi (newval); 
do 
-{ 

if(*newval  <  'O'  |!  ‘newval  > 
if(*newval  !*  '.') 


fprintf(stderr,  "Invalid  input  -  enter  nunerical  value  \n"); 
retum(-l) ; 


— }  v^iile(*(++newval)  !*  '\0'] 

if(num  <011  nun  >  OxFFFF) 
retum(-l) ; 
return (num) ; 
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************************************************************************ 


/*  g  fnum.c  */ 
/*  g_fnum.c  Prompt  user  for  floating  point  numerical  input  and  */ 
/*  then  verify  it.  If  valid  string  is  found,  returns  length  of  number  */ 
/*  string,  otherwise  returns  zero.  */ 
/*  created  6/6/86  B.  Lentz  */ 
/*  */ 


j *************************************************** ********************* j 

♦include  <stdio.h> 

double  g  fnum(newval) 
char  ‘newval; 

r{ 

int  i;  /*  counter  for  loop  */ 

float  fnum;  /*  nunerical  value  of  parameter  */ 

double  atof(); 

char  string[20];  /*  input  string  containing  new  value  */ 

printf (“\n\tEnter  new  value:  "); 
i f (fgets (string,  20,  stdin)  *•  0) 
printf ("\nerror  in  input\n"); 
sscanf(string,  "%s",  newval); 
fnum  »  atof (newval); 
do 

if(*newval  <  ‘O’  1!  ‘newval  >  '9') 
if (‘newval  !* 

— { 

fprintf(stderr,  "invalid  input  -  enter  numerical  value  \n"); 
retum(-l); 

L-} 

— }  while(*(++newval)  !*  ' \0 ' ) ; 
return  (f  nun) ; 
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j  I  ft***********************************************************************  j 

2  /*  g  parm.c  */ 

3  /*  g_parm.c  reads  text  file  containing  parameters  into  ’usr1'  structure  */ 

4  /*  created  4/1/85  B.  Lentz  */ 

5  /*  */ 

6  /*  modification  3/13/86  added  date  parameter  */ 

7  /*  modified  6/12/86  rt's  */ 

8  /*  V 

g  /  ********ih*  ★**★*★*★★★*********★**★*★*★★*★**★*******★***★*★*  **♦**★**■**★*★*  j 

10  #include  <stdio.h> 

11  linclude  "parm.h" 


13  g_parm(usr,  s3) 

14  struct  value  *usr; 

15  struct  value  *s3; 


/*  user  value  */ 
/*  S-CUBED  */ 


16 

*7  r-{ 


18  FILE  *pfile;  /*  file  pointer  to  text  file  containing  parameters  */ 

19  int  j;  /*  coulter  */ 

20  char  keyword [16]; 

21  char  string[20];  /*  temporary  string  storage  for  comparison  */ 

22  char  str2[20];  /*  temporary  string  storage  for  comparison  */ 

23  int  badparm;  /*  flag  for  valid  parameter  */ 


22  char  str2[20]; 

23  int  badparm; 


25  if  ((pfile  *  fopen(targfi le[0] ,  "r"))  **  0) 


fprintf (stderr,  "\nERR0R  opening  parameters  file  %s\n",  targf i le[0] ) ; 
exit(0); 


30  while((fscanf (pfile,  ”%s",  keyword))  !*  0) 


if(strcmp(keyword,  "runtime")  **  0) 

1  f$canf(pfile,  "%f  %f",  &usr->runtime,  Ss3->runtime); 


else  if (strcmp(keyword,  "caldate")  **  0) 
r~{ 

fscanf (pfile,  "%s  %s",  &usr->caldate,  &s3-xaldate); 

M 

else  if (strcmp(keyword,  "date")  **  0) 


fscanf (pfile,  "%s  %s",  &usr->date,  &s3->date) ; 
— } 

else  if(stronp(keyword,  "ca-^ag")  ■*  0) 


L  fscanf (pfile,  "%d  %d" ,  &usr->calflag,  &s3->calflag); 

} 

else  if (strcmp(keyword,  "spikcnt")  **  0) 

H 

|  fscanf(pfi le,  "%d  %d",  Susr->spikcnt,  &s3->spikcnt) ; 
1— l 


else  if(strcmp(keyword,  "verbose”)  **  0) 


fscanf (pfile,  "%s  %s",  keyword,  string); 
usr->verbose  *  logstr(keyword) ; 


R 
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$3->verbose  1  logstr(string) ; 


-1 


else  if  (strcmp(keyvord,  "gain')  ==  0) 

-{ 

fscanf(pfi le,  "%u  %u",  &usr->gain,  8s3->gain); 
-} 

else  if (strcmp(keyword,  "freq")  *»  0) 
fscanf(pfi le,  "%f  %f",  4usr->freq,  8s3->freq); 


-} 


0) 


else  if (strcmp( keyword,  "smooth") 
r-{ 

fscanf(pfile,  "%d  %d",  8usr->smooth,  8s3->smooth) ; 

L-} 

else  if (strcmp(keyword,  "threshl")  —  0) 


-{ 


fscanf(pfile,  "%d  %d",  8usr->thresh[l] ,  8s3->thresh[l]); 


else  if (strcmp(keyword,  "thresh2") 


0) 


-{ 


fscanf{pfile,  "%d  %d",  8u$r->thresh[2] ,  8s3->thresh[2]) ; 


— } 


else  if (strcmp(keyword,  "thresh3")  **  0) 


-{ 


fscanf(pfile,  "%d  %d",  &usr->thresh[3] ,  8s3->thresh[3]); 


0) 


else  if (stremp( keyword,  "thresh4") 

-{ 

fscanf(pfile,  "%d%d",  8usr->thresh[4] ,  8s3->thresb[4j) ; 


-} 


0) 


else  if (strcmp( keyword,  "minarea"' 
r~{ 

fscanf(pfile,  "%d  %d",  &usr->ninarea,  &$3->mi narea); 
-} 

else  if (strcmp (keyword,  "nointstd")  **  0) 

-{ 

fscanf(pfi le,  "%d  %d",  8usr-*istd,  8s3->nstd); 

1 — } 

else  if (strcmp(keyword,  "nununks")  »»  0) 

-{ 

fscanffpfile,  "%d  %d",  8usr->nunks,  &s3-^nunk5); 


else  if (strxmp( keyword,  "nlevels") 


0) 


fscanf (pf i le,  "% d  %d",  8usr->n levels,  8$3->n1evels); 

-} 

else 


badparm  «  TRUE; 

for(j  ■  1;  j  <*  &usr->nstd;  j++) 
-{ 

sprintf (string,  ’stdRT%d' ,  j); 
sprintf (str2,  "stuarea%d",  j); 
if (Stroup (keyword,  string)  **  0) 

H 
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fscanf(pfile,  '  %d  %d" ,  &usr->stdrt[j] ,  &s3->stdrt[j]); 
badparm  *  FALSE; 
break ; 

else  if(strciH3(keyword,  str2)  -«  0) 

p{ 

fscanf(pfi le,  "%d  %d",  &usr->stdarea[j],  8s3->stdarea[j]) ; 

badparm  •  FALSE; 

break; 

L-} 

— }  /*  end  for  j  */ 
if (badparm) 

fprintf(stderr,  "\nParameter  %s  is  not  used\n",  keyword); 

— }  /*  end  else  */ 

— }  /*  end  while  */ 
fclose(pfile) ; 

usr->period  »  1000000/ (usr->freq  *  2.5  *  8.0); 
usr->mjnpt  *  60  *  usr->freq  *  usr->runtime; 
usr->nconv  *  usr->nurnpt  *  8; 
for  (j  *  1;  j  <■  usr->n$td;  j++) 

r{ 

usr->sumarea[j]  *  0; 
usr->sunrt[j]  *  0; 

M 

— }  /*  end  g  struct  */ 


logstr(string) 
char  ‘string; 

p{ 

if(strcmp(string,  "ON")  *«  0) 
return (TRUE) ; 
else 

return (FALSE) ; 


09-11-87  08:40:20  g_sspec.c 
Fri  09-11-87  10:04:51  g  sspec 


pg 

of  1 
1-55 


y *********************************************** ************ ********* ****y 
/*  g  sspec. c  */ 

/*  g  sspec. c  reads  text  file  containing  target  parameters  into  V 

/*  structure  "samp"  */ 

/*  created  4/1/85  B.  Lentz  */ 

/*  */ 

y************************************************************************ y 

linclude  <stdio.h> 
linclude  "parm.h" 
linclude  "targ.h" 


g_sspec(samp) 
struct  sample  sampQ ; 


FILE  ‘pfile;  /*  file  pointer  to  text  file  containing  parameters  */ 

int  i,  j;  /*  coulter  */ 

int  badparm;  /*  flag  incorrect  parameter  */ 

char  keyword [16] ;  /*  name  from  text  file  */ 

char  string[20];  /*  holds  parameter  name  */ 

extern  struct  value  usr; 


for(i  »  1;  i  <■  usr.nunks;  i++) 


if  ((pfile  *  fopen(targfile[i],  " r" ) }  ■*  0) 

-{ 

fprintf(stderr,  "\nERR0R  opening  parameters  file  %s\n",  targfile[i]); 
exit(0) ; 

-} 

while((fscanf (pfile,  "%s",  keyword))  !*  0) 

-{ 

if(strcmp(keyword,  "rettime")  0) 

[— < 

fscanf(pfile,  "%d  %*d",  &samp[i].rt); 

— } 

else  if(strcmp(keyword,  "name")  **  0) 
fscanf (pfile,  "%s  %‘s",  &samp[i].name); 

L-} 

else  if (stranp(keyword,  'units")  ■*  0) 


fscanf (pfile,  "% s  %*s",  &samp[i], units); 


else  if (strcmp (keyword,  "spike")  ■*  0) 
j  fscanf (pfile,  "%f  %*f",  Ssamp[i].  spike); 
else  if ( strcmp (keyword,  'precisn")  0) 


fscanf(pfi!e,  "Vf  %*f",  4samp[>]. precisn); 


else  if (strcmp (keyword,  accurcy')  **  0) 
fscanf(pfile,  ' %f  %*f",  Ssampfi]. accurcy); 


*  xV/.  V 


s  s 
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else  if(strcmp(keyword,  "cutoff)  **  0) 


-{ 


| _ fscanf (pf ile,  "%f  %*f,  &samp[i]  .cutoff)  ; 

else  if(strcmp(keyword,  "qtlimit")  •*  0) 


-{ 


j  fscanf(pfile,  "%f  %*f",  &samp[i] -qtlimit) r 
* — } 

else  if(stranp(keyword,  "dilution")  **  0) 


-{ 


Lfscanf  (pfile,  "%f  Vi",  &samp[i]. dilution); 

} 

else 
i — { 

badparm  *  TRUE; 

for(j  *  1;  j  <*  usr.nstd;  j++) 
r~{ 

sprintf(string,  "stdconcW,  j); 
if(strcnp(keywond,  string)  **  0) 

fscanf(pfi1e,  "%f  V f",  &samp[i].conc[j]); 
badparm  *  FALSE; 
break ; 

-} 

sprintf (string,  "slope%d",  j); 
if(strcmp(keyword,  string)  **  0) 

fscanf(pfile,  "%f  %*f",  &samp[f] .slopefj]) ; 
badparm  =>  FALSE; 
break ; 

-} 

sprintf(string,  "intercp%d",  j); 
if(strcmp(keywond,  string)  *»  0) 

'  f$canf(pfi le,  “Hi  %*f”,  &sanp(i].intercp[j]); 
badparm  *  FALSE; 
break; 

L_} 

-} 

if(badparm) 

fprintf (stderr,  "parameter  %$  is  not  used\n".  keword) ; 

M 

L-) 

if (fclose(pf i le)  **  -1) 

fprintf (stderr,  "\nERR0R  closing  %s\n",  targf ile[i]) ; 
sampfi]  .sumurt  *  0; 
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j *******************th>*************** ***************************** *******/ 

/*  initdio.c  */ 

/*  set  digital  I/O  bit  0  of  port  0  for  output  and  initialize  to  open  */ 
/*  relay  switch  (set  to  1)  */ 

/*  this  state  will  cause  the  Tekmar  to  wait  for  the  switch  closure  to  */ 
/*  begin  desorb  */ 

/*  created  2/7/86  B.  Lentz  */ 

/*  */ 

/ft*********************************************************************** / 


linclude  <$tdio.h> 

linclude  "bitset.h"  /*  bit  setting  macros  for  data  acquisition  board  */ 

initdioQ 

-{ 

unsigned  char  status;  /*  used  for  error  check  */ 
if  ((STAT  REG  &  0x70)  !*  0) 


fprintf(stderr,  "\nFATAL  ERROR-1 1  legal  status  register  va1ue\n"); 
fprintf (stderr,  "\nStatus  Register  value  is  sso\n" ,  STAT_REG); 
exit(O) ; 

} 


COMMREG(CSTOP) ; 
status  ■  0ATA_0UT; 
while(!(STAT_REG  &  COMM  WAIT)); 
COMMREG(CCLEAR); 

whi  le( !  (STAT_REG  &  CCTM  WAIT)) ; 
COMMREG(CSOUT); 

whi le(STAT_REG  &  WRITE  WAIT); 
0ATA_IN(0I0P0RT) ; 

while(!(STAT_REG  &  C0MM_WAIT)); 
status  ■  STAT_REG; 
if (status  &  0x80) 
ioerr(); 

setdio(0,  CDIOOUT) ; 

-} 


I 
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j  j **************** ******************************************************** i 


2  /*  ioerr.c  */ 

3  /*  read  status  registers  for  error  reporting  from  data  translation  */ 

4  /*  board  */ 

5  /*  created  2/6/86  B. Lentz  */ 

6  /*  */ 


^  /*******************  ***************  ************************  **************  j 

8 

9  #include  <stdio.h> 

10  #include  "bitset.h" 

11 

12  ioerr() 

13  r{ 

14  unsigned  cbar  temp; 

15  char  errorl; 

16  char  error2; 

17 

18  fprintf(stderr,  "\nFATAL  BOARD  ERROR  \n"); 

19  fprintf(stderr,  "\nStatus  Register  Value  is  %x\n",  STAT_RE6) ; 

20  /*  read  error  register  */ 

21  COmREG(CSTOP)  ; 

22  temp  -  DATAJXJT; 

23  whi  le(!  (STAT_REG  &  C0MM  WA1T) ) ; 

24  C0MM_REG(CERR0R) ; 

25  whi le( ! (STAT_REG  &  READ_WA1T) ) ; 

26  errorl  -  DATAJXJT; 

27  whi  le( !  (STAT_REG  &  READ  WAIT) ) ; 

28  error2  *  DATAJXJT; 

29  fprintf (stderr,  "Error  Register  values  are:  \n"); 

30  fprintf(stderr,  "bytel:  %x  \tbyte2:  %x",  errorl,  error2); 

31  exit(O) ; 

32  I—} 

33 


/*  first  byte  of  error  code  */ 
/*  second  byte  of  error  code  */ 
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/*  match. c 

/*  routine  to  locate  standard  peaks  (internal  and  target) 

/*  within  standard  samples  and  determine  actual  retention  times  and 
/*  response  factors  for  each  standard. 

/*  created  7/11/85  B.  Lentz 


♦include  "panm.h" 
♦include  "targ.h" 
♦include  <stdio.h> 

float  match(npeaks, 
int  npeaks; 
int  rawrtQ ; 
float  rwaneaQ; 
long  finareaO  ; 
int  rtimeQ; 


raw rt,  rwarea,  finarea,  rtime) 

/*  total  number  of  peaks  found  */ 

/*  retention  times  of  all  peaks  */ 

/*  areas  of  all  peaks  */ 

/*  areas  of  selected  peaks  */ 

/*  retention  times  of  selected  peaks  */ 


extern  struct  sample  samp[] 
extern  struct  value  usr; 
int  windowl;  /* 

int  wind<*£;  /* 

int  start [2];  /* 

int  end[2];  /* 

int  index;  /* 

int  i,  j,  k,  m:  /* 

int  rrt;  /* 

float  factor;  /* 

float  nun;  /* 

double  sqrt();  /* 

float  fstarea;  /* 

float  farea;  /* 

float  den;  /* 

float  uden;  /* 

float  rden;  /* 

float  match;  /* 

float  hit;  /* 

float  totpts;  /* 

int  fouidpk;  /* 

int  tenp[MAXSTDS] ;  /* 


;  /*  target  cmp  parameters  */ 

/*  contains  parameters  */ 
range  in  which  to  search  for  std  peak  */ 
range  in  which  to  search  for  std  peak  */ 
beginning  of  search  range  V 
end  of  search  range  */ 
index  to  start  of  target  parameters  */ 
loop  counters  */ 

relative  retention  times  based  on  first  std  */ 

ratio  of  actual  retention  time  to  estimate  */ 

numerator  of  match  factor  */ 

square  root  function  V 

std  area  expressed  as  floating  pt/1000  V 

actual  area  expressed  as  floating  pt/1000  */ 

denominator  of  match  factor  */ 

portion  of  den  containing  sample  contributions  */ 

reference  denominator  */ 

factor  indicating  quality  of  match  V 

best  match  factor  found  V 

total  number  of  data  points  collected  */ 

boolean  used  to  flag  potential  match  V 

temporary  storage  for  index  to  peak  */ 


hit  -  0.0; 

totpts  •  usr. runtime  *  usr.freq  *  60; 

window2  *  (int)  (totpts  *  0.03);  /*  3%  window  */ 

windowl  •  (int)  (totpts  *  0.40);  /*  40%  window  */ 

start [0]  ■  usr.stdrt(l]  -  windowl; 
if(start[0]  <  0) 
start[0]  *  0; 

end[0]  •  usr.stdrtfl]  +  windowl; 
for(i  *  1;  i  <*  npeaks;  i*+) 

i — { 


^  v  *t» «r. tv V. .'  .v '„■. .V.'-V- '.V  V.  V V 


09-11-87  08:41:58  match. c 
Fri  09-11-87  10:04:51 


match 


>v. 


-\  - 
-\‘- 

fl- 

*;■ 


v- 

r.j 


I 


56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
73 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 


if((rawrt[i]  >  start [0])  &&  (rawrt[i]  <  end[0])) 


-{ 


if(rwarea[i]  >  usr.minarea) 

•{ 

temp[l]  *  i; 

factor*  (float)  rawrt[i]/(usr.stdrt[l]); 

for(k*2;  k  <*  usr.nstd;  k-n-) 

I — { 

temp[k]  *  0; 

rrt  *  (int)  (usr.stdrt[k]  *  factor); 
start[l]  *  rrt  -  windix2; 
i f (start [1]  <  0) 
start [1]  *  0; 
end[l]  *  rrt  +  window2; 

for(j  *  1;  j  <*  npeaks;  j++) 

'  if ((rawrt [j]  >  start[l])  &&  (rawrt[j]  <  end[l))) 
-{ 


if(rwarea[j]  >  usr.minarea) 
temp[k]  *  j; 


found pk  *  TRUE; 

for(m  -  1;  m  <*  usr.nstd;  im+) 

if (tanp[m]  »*  0) 
foundpk  *  FALSE; 

-} 

if (foundpk) 

r~{ 

nun  *  0.0; 
uden  *  0.0; 
rden  *  0.0; 

for(m  *  1;  m  <*  usr.nstd;  nr+) 

I — { 

fstarea  *  (usr.stdarea[m])/1000. ; 
farea  *  rwarea[tOTp[m]]/1000. ; 
num  +*  (fstarea  *  farea); 
uden  +*  (f area  *  f area); 
rden  +*  (fstarea  *  fstarea) ; 

' — } 

if  ((den  *  sqrt(uden  *  rden))  **  0.0) 

l — (  . 

printf ("\n\n  ERROR  -  denominator  is  zero\n  ); 

match  *  0.0; 

M 

else 

r- ^ 

match  *  num/den; 

' — } 

i f (match  >  hit) 
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/*  menutest.c  User  menu  */ 

/*  This  is  the  main  menu  routine  to  control  sample  automation  on  the  gc  */ 
/*  Created  3/28/85  last  revision  4/21/86  */ 

/*  B.  Lentz  */ 

/A***********************************************************************/ 

linclude  <stdio.h> 
fdefine  TRUE  (0--0) 

Idefine  FALSE  ITRUE 


main(argc,  argv) 
int  argc; 
char  *argv[] ; 

-{ 

int  instmo; 


/*  menu  selection  */ 


while(0**0)  /*  infinite  loop  */ 

r-{ 

system("cls '") ;  /*  clear  screen  */ 

printf ("\n\n\n\n\n\n\t\tSe1ect  one  of  the  following  options:\n“); 
printf("\n\n\t\t  1)  Analyze  for  TCE  \n"); 
printf ("\t\t  2)  Calibration  run  for  TCE\n"); 
printf ("\t\t  3)  Prepare  new  Data  diskette(for  drive  B)\n"); 
printf ("\t\t  4)  Prepare  new  V0A  Data  System  diskette(for  drive  A)\n"); 
printf ("\t\t  5)  Retrieve  QC  data\n"); 
printf ("\t\t  6)  Retrieve  archived  data\n"); 
printf ("\t\t  7)  Change  operating  parameters^” ) ; 
printf ("\t\t  9)  Exit  Menu\n"); 
printf ("\n\n\t  ENTER  OPTION  NLWBER:\t") ; 

/*  get  input  from  keyboard  */ 
instmo  *  selectQ  ; 
whi  le( instmo  <  1  [|  instmo  >  9) 

r{ 

fprintf (stderr,  "\nlnvalid  Input.  Enter  NUMBER  of  selected  option:\t"); 
instmo  *  select (); 

M 

system("cls") ; 
switch(instmo) 

r{ 

case  1:  /*  sample  run  */ 

r{ 

printf ("Analyze  for  TCE\n" ) ; 
system("analyzl  0") ; 
break ; 

— }  /*  end  case  2  */ 
case  2:  /*  calibration  run  */ 

r-{ 

printf ("Calibration  run  for  TCE\n"); 
system("analyzl  1"); 
break ; 

L> 

case  3:  /*  format  data  diskette  *! 


|  printf(" 


Prepare  data  diskette  for  use\n"); 


I. 
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system("fmtdis ") ; 
break ; 

M 

case  4:  /*  Format  system  floppy  disk  */ 

r{ 

printf("\nPrepare  new  VOA  Data  System  diskette\n") ; 

system("copsys") ; 

break; 

L} 

case  5:  /*  QC  */ 

r-{ 

printf ("Retrieve  QC  Data\n"); 

system("qcexec ") ; 

break; 

«-} 

case  6:  /*  retrieve  data  from  floppy  */ 

r{ 

printf ("Retrieve  archived  data  from  diskette\n"); 

system("xretrv") ; 

break; 

case  7:  /*  Modify  operating  parameters  */ 

r{ 

printf ("Change  operating  parameters'^ ") ; 

system("change") ; 

break; 

L} 

case  9:  /*  exit  to  DOS  */ 

r{ 

printf ("\nReturning  to  operating  system\n"); 

printf ("Type  \"demo\"  to  return  to  data  analysis  system\n"); 

exit(0); 

break; 

*-} 

default: 

r-{ 

printf ("\nNot  a  valid  selection\n"); 
break; 

-} 


■r.V-V'y. 
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/*★★**★*****★* ****************************************************** *****/ 


/*  pkstart.c  */ 
/*  Routine  called  by  pk  det.c  to  determine  it  a  new  peak  */ 
/*  has  been  found.  */ 
/*  created  5/14/85  B.  Lentz  V 
/*  V 


I *****************************************♦*****»************************/ 

♦define  TRUE  (0  =*  0) 

♦define  FALSE  ! TRUE 


pkstart(np,  thrl,  thr2) 
int  npQ  ; 
int  thrl; 
int  thr2; 

r{ 

int  i ; 

if((np[l]  -  2*np[2]  +  np[3])  >0) 
r-{ 

if((np[l]  -  3*np[3]  +  2*np[4])  >0) 
if((np[l]  -  4*np[4]  +  3*np[5])  >0) 

/*  check  for  false  peak  start  */ 
if ( (np[5]  -  np [1] )  >  thrl) 

r{ 

for(i  »  1;  i  <  4;  i++) 

if(((np[i+2]  -  np[i+l])  -  (np[i+l]  -  np[i]))  >*  thr2) 
return (TRUE); 

*-} 

if(((np[5]  -  np[4] )  -  (np[2]  -  np[l] ) )  >•  thr2) 
retum(TRUE) ; 
else 

retum(  FALSE) ; 

— }  /*end  false  peak  start  */ 
return(FALSE)  ; 

I—} 
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/*  pk_det.c  *! 
/*  Main  section  of  peak  detection  software.  Initializes  data  */ 
/*  segments  and  calls  subroutines  to  establish  peak  start,  peak  end  */ 
/*  and  peak  max.  *! 
/*  created  5/14/85  8.  Lentz  */ 
/*  V 
!*+** ************************************************************* **+*++*/ 
♦include  <stdio.h> 

♦include  "parm.h" 


♦define  DERI VI (a,  b,  c) 
♦define  DERIV2(a,  b,  c) 

static  int  x$t [NPEAKS] ; 
static  int  yst [NPEAKS] ; 
static  int  xend[NPEAKSj; 
static  int  yend [NPEAKS] ; 
static  float  baseb; 
static  float  basera; 


((-3*a  +  4*b  • 
(a  -  2*b  +  c) 


0/2) 


/*  starting  x  value  for  peak  */ 

/*  starting  y  value  for  peak  */ 

/*  value  of  x  at  peak  end  V 
/*  value  of  y  at  peak  end  */ 

/*  intercept  of  baseline  segment 
/*  slope  of  baseline  segment  */ 


pk  det(npts,  data,  usr,  xmax,  area) 


unsigned  npts; 
unsigned  int  ‘data; 
int  xmax [NPEAKS] ; 
float  areafNPEAKS] ; 
struct  value  *usr; 


/*  nunber  of  data  points  */ 

/*  data  array  */ 

/*  value  of  x  at  peak  max  */ 

/*  total  peak  area  */ 

/*  paramters  including  threshold  values 


H 


’nt  l,  j; 
int  np[8] ; 
int  pos; 
int  status; 
int  ct; 
irit  peak; 
int  pkstart(); 
int  integrQ; 
int  pkcnt; 


/ 


/*  subset  of  data  V 
/*  current  x- value  of  data  */ 

/*  was  conditional  satisfied?  */ 

/*  count  of  times  condition  is  met 
/*  was  peak  identified?  */ 

/*  subroutine  to  locate  start  of  peak 
/*  integrates  peak  areas  */ 

/*  nuiber  of  peaks  detected  */ 


int  ymax [NPEAKS] ;  I*  value  of  y  at  peak  max  */ 


int  rd_data(); 
int  pk_end(); 
float  TT ; 
int  end; 
int  baseline; 
long  base; 
int  bent; 
int  bm; 
int  mp; 

pkcnt  *  0; 
npts  —  5; 


/*  subroutine  to  read  next  segment  of  data  */ 


/* 

/* 

/* 

/* 

/* 

/* 


subroutine  to  predict  peak  end  */ 
distance  to  predicted  peak  end  */ 
reached  peak  end?  */ 
current  value  of  baseline  V 
sun  baseline  points  */ 

♦  of  pts  used  to  determine  baseline  ‘ 
distance  from  peak  start  to  peak  max 


/ 

V 


/* 

/*  distance  from  peak  max  to  predicted  oeak  end  */ 

I*  nunber  of  peaks  in  cluster  */ 

/*  examine  5  point  segments  */ 


/*  initialize  data  segment 
for  { i »2 ;  i<>;  i++) 


V 


r' 


np[i]  *  data[ i -2] ; 
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xst[i]  *  0; 
yst [i ]  *  0; 
XDHx[i]  *  0 
ymax[i]  •  0 
xend[i]  *  0 
yend[i]  •  0 
area[i]  ■  0.0; 
-} 

posO; 

j-4; 

while(pos<iipts) 


base  *  0; 
bent  -  0; 

whi1e(!(peak  *  pkstart(np,  usr->thresh[l] ,  usr->thresh[2]))) 

{ 

pos  *  rd  data(np,  pos,  data); 
if(pos  >  npts) 

i — { 

if  (pkent  <  1) 
retum(pkcnt); 
integr(pkcnt,  area,  data); 
if(usr->verbose) 

pk_prt(pkcnt,  xmax,  ymax,  area); 
return (pk ent); 

-> 

bcnt++; 

base  +*  np[l]; 
baseline  *  base/bent; 

-} 

newpk:  pkcnt++; 
if (pkent  >  NPEAKS) 

pkent—; 

printf("\n  exceeded  max  nutiber  of  peaks\n"); 
integr(pkcnt,  area,  data); 
if (usr->verbose) 

pk  prt(pkcnt,  xmax,  ymax,  area); 
return(pkcnt) ; 

area[pkcnt]  *  (float)  np[lj; 
yst  [pkent]  *  np[l]; 
ymax[pkcnt)  *  np[l]; 
xmax [pkent]  *  pos; 
xst [pkent]  *  pos; 

ct  *  1; 

status  *  FALSE; 

I — d°  ( 

if(0ERlVl(np[l],  np[2],  np[3])  <  0) 
if  (DERIV2 (np[l] ,  np[2] ,  np[3])  <  1) 

pi 

|  if(usr->verbose) 

48 
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printf("shoulder  at  pos  HA  \n",  pos); 
status  *  TRUE; 
if (ct  >  1) 

-{ 

ct  *  pos  -  ct/2; 
xmax [pkcnt]  •  ct; 
ymax[pkcnt]  ■  data[ct] ; 

-} 

} 


pos  *  rd  data(np,  pos,  data); 
if (pos  >  npts) 

pkcnt--; 

integr(pkcnt,  area,  data); 
if(usr->verbose) 

pk_prt( pkcnt,  xmax,  ymax,  area); 
retum(pkcnt); 

*-} 

if(np[l]  —  ymax[pkcnt])  /*  flat  area  */ 

ct++; 

else  if  (np[l]  >  ymax  [pkcnt]) 

r~{ 
ct  •  1; 

ymax  [pkcnt]  -  np[l]  ; 
xmax [pkcnt]  »  pos; 

if((np[2]  <  np[l] )  &&  (np[3]  <  np[2])  &&  (np[4]  <  np[3])  &&  (np[5]  <  np[4])) 
r-{ 

status  •  TRUE; 
if (ct  >  1) 
r-{ 

ct  -  pos  -  ct/2; 
xmax [pkcnt]  -  ct; 
ymax  [pkcnt]  ■  data[ct]; 

L_} 

*— } 

area[pkcnt]  +■  (float)  np[l] ; 

— }  while  ( ! status) ; 

bm  *  xmax[pkcnt]  -  xstfpkcnt] ;  /*  peak  start  to  apex  V 

— do  { 

if(DERIVl(np[l] ,  np[2],  np[3])  >0) 

— { 

if (0ERIV2 (np[l] ,  np[2],  np[3])  >0) 
r-{ 

if(usr->vertose) 

printf("shoulder  at  pos  %d  \n",  pos); 
xend [pkcnt]  *  pos; 
yend [pkcnt]  »  np[|] ; 
goto  newpk; 

*-} 

pos  *  rd_data(np,  pos,  data); 
if (pos  >  npts) 


r *>  V •>  V  Vjr>> v V V 
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pkcnt--; 

integr(pkcnt,  area,  data); 
if(usr->verbose) 

pk_prt(pkcnt,  xmax,  ymax,  area); 
retum(pkcnt); 

*-} 

area[pkcnt]  +■  (float)  no[l]; 

if(peak  *  pkstart(np,  usr->thresh[l] ,  usr->thresh[2])) 
i — { 

xend [pkcnt]  *  pos ; 
yend [pkcnt j  *  n p[l]; 
goto  newpk; 

I — } 

if((np[4]  -  np[5])  <  usr->thresh[4]) 
end  ■  TRUE; 
else 

end  *  pk_end(np,  usr->thresh[3]); 
if((np[5]  -  baseline)  >  (ymax[pkcnt)  -  baseline) *0.75) 
end  ■  FALSE; 
if(np[5]  <  baseline) 
end  •  TRUE; 

mp  *  (pos+4)  -  xmax[pkcnt]; 
if(mp  >  2*bm) 
end  «  TRUE; 

— }  while('.end); 

if(mp*2  <  bm) 

H  -  bm*2; 
else 

TT  *  mp*2; 
end  *  FALSE; 
ct  *  0; 
do 
I — { 

if(DERIVl(np[l] ,  np[2],  np[3])  >0) 
i — i 

if (DERIV2(np[l] ,  np[2],  np[3])  >0) 
r—{ 

if(usr->verbose) 

printf("shoulder  at  pos  %d\n",  pos  ); 
xend [pkcnt]  ■  pos; 
yend  [pkcnt]  •  np(l]; 
goto  nev^jk; 

L-} 

I — } 

pos  *  rd  data(np,  pos,  data); 
if (pos  >  npts) 
i — { 

pkcnt—; 

integr(pkcnt,  area,  data); 
if(usr->verbose) 

pk  prt (pkcnt,  xmax,  ymax,  area); 
return (pkcnt); 
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pk_prt(pkcnt, 
int  pkcnt; 
int  xmaxQ; 
int  ymax[] ; 
float  areaQ ; 


xmax,  ymax,  area) 

/*  number  of  peaks  detected  */ 
/*  value  of  x  at  peak  max  */ 

/*  value  of  y  at  peak  max  */ 

/*  total  peak  area  */ 


I  int  k; 


printf("\n\t  start\t\t  center\t\t  end  \t\t  area  \n"); 
for(k*l;  k  <■  pkcnt;  k++) 

— { 

printf("  x  %d  y  %d\t",  xst[kj,  yst[k]); 
printf("  x  %d  y  %d\t",  xmax[k],  ymax[k]); 
printf("  x  %d  y  %d\t",  xend[k],  yend[k]); 
printf("  %f\n" ,  area[k]); 

L-} 

l—} 


/*  integr.c  used  to  integrate  peak  areas  */ 
/*  created  5/28/85  */ 
/*  calls:  segm  */ 
/*  B.  Lentz  */ 


/ ♦+**  +  *•*********«**********+****★***********«******«********★*****★*******★*  j 


integr(pkcnt,  area,  data) 

int  pkcnt;  /*  nunber  of  peaks  detected  */ 

float  area[] ;  /*  total  peak  area  */ 

unsigned  int  ‘data; 


int  start; 
int  end; 
int  y; 
int  pos; 
int  j ; 


/*  array  position  in  peak  cluster  */ 

/*  array  position  at  end  of  cluster  */ 

/*  y  value  from  baseline  */ 

/*  current  array  position  in  peak  cluster  */ 


long  y 1 ,  y2 ;  /*  baseline  values  at  start  and  end  of  peak  */ 

long  sun;  /*  area  below  peak  baseline  V 

int  segm(); 

extern  struct  value  usr; 


start  «  1; 
end  «  pkcnt; 
while(start  <»  pkcnt) 

— { 

end  *  segm(start,  end,  data); 
if(usr. verbose) 

printf ( "\tsegst  %d  segend  *d\n",  xst[start],  xend[end]); 
for(j  *  start;  j  <•  end;  j-h.) 
r~{ 

if((yl  •  (xst[j]  *  basem)  ♦  baseb)  *•  0) 
yl  *  0; 

if((y2  *  (xend[j]  *  basem)  ♦  baseb)  •«  0) 

52 


Pg  51 
of  105 
276-330 


->CV  V*  TV  W  V  ,’VCTA.V'  ^JV/W^JWTV'.V'JVTUL’B 


09-11-87  08:45:00  pk  det.c 
Fri  09-11-87  10:04:51 


Pg  52 
of  105 
331-385 


y 2  ■  0; 

sum  *  (xend[j]  -  xst[j]  +  1)  *  ((yl  +  y2)/2); 
area[j]  -«  (float)  sum; 

-} 

start  »  end  +  1; 
end  •  pkcnt; 

} 


ft**************************************************************************/ 

/*  segm.c  recursive  routine  used  to  set  local  baselines  */ 

/*  created  6/12/85  V 

/*  B.  Lentz  */ 

Z^*******^********^*)******^*!^*****************************************^ 


segm(start,  end,  data) 

int  start;  /*  array  element  to  start  baseline  */ 

int  end;  /*  array  element  to  end  baseline  V 

unsigned  int  ‘data; 


float  mck; 
float  fxl; 
float  fyl; 
float  fx2; 
float  fy2; 
int  a; 
int  yl,  y2; 
int  x; 


/*  check  slope  value  */ 


/*  current  peak  position  »/ 

/*  values  from  baseline  */ 

/*  current  baseline  position  */ 


fxl  ■  xst [start]; 
fx2  ■  xend[end]; 
fyl  ■  yst  [start]; 
fy2  •  yend[end] ; 


baseb  •  ((fx2  *  fyl)  -  (fxl  *  fy2))/(fx2  -  fxl); 
basem  •  (fyl  -  baseb)/fxl; 
mck  *  (fy2  -  baseb)/fx2; 
for(a  *  end;  a  >•  start;  a—) 


yl  •  (xst [a]  *  basem)  +  baseb; 
y2  •  (xendfa]  *  base*)  +  baseb; 
if  ((yl  >  yst[a] )  ||  (y2  >  yend[a])) 


a—; 

if (a  start) 

H 

fx2  *  xend[start]; 
fy2  ■  yend [start]  ; 

baseb  ■  ((fx2  *  fyl)  -  (fxl  *  fy2))/(fx2  -  fxl); 
basem  ■  (fyl  -  baseb) /fxl; 
return (start); 

-} 

end  *  segm($tart,  a,  data); 
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else  i f (a  >  start) 

if(xst[a]  >  xend[a-l])  /*  baseline  between  peaks  */ 
for(x  *  xend[a-l];  x  <  xst[a];  x++) 

r{ 

/*  check  for  points  below  chromatographic  baseline  */ 
yl  *  (x*basem)  +  baseb; 
if (yl  >  data[x]) 


a—; 

if (a  <*  start) 

-{ 

/*  segment  only  includes  one  peak  */ 
fx2  *  xend [start] ; 
fy2  *  yend  [start]; 

baseb  -  ((fx2  *  fyl)  -  (fxl  *  fy2))/(fx2  -  fxl); 
basem  »  (fyl  -  baseb)/fxl; 
retum(start)  ; 

-} 

end  *  segm(start,  a,  data); 

-} 


return (end) ; 
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/*************************<M********<H*************************************  j 

!*  pk_end.c  */ 
/*  Predict  peak  end  for  pk  det.c  */ 
/*  created  5/15/85  B.  Lentz  V 


/********************************************************★*************** / 


♦define  TRUE  (0  0) 

♦define  FALSE  !TRU£ 

pk_end(np,  thr) 
int  npQ; 
int  thr; 


1;  i  <  4;  i++) 


if(((np[i+2]  -  np[i+l])  -  (np[i+l]  -  np[i]))  >-  thr) 
retum(  FALSE); 


if(((np[5]  -  np[4] )  -  (np[2] 
return  (FALSE)  ; 
else 

retum(TRUE)  ; 


np[l]))  >-  thr) 
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j a********************************* **** ******* *************************** 


/*  pmparm.c  */ 
/*  Print  current  user  defined  parameters  */ 
/*  created  8/14/85  B.  Lentz  */ 

/*  V 


j **************************************************** ********************** j 

linclude  "parm.h1 

pmparm(usr) 
struct  value  *usr; 

|-{ 

int  i; 

float  convrt;  /*  convert  from  array  position  to  minutes  */ 

system(  els") ; 
convrt  *  usr->freq  *  60.0; 
if(usr->verbose  --  TRUE) 
printf("\n\n\na)  Verbose  option  is  0N\n' ) ; 
else 

printf("\n\n\na)  Verbose  option  is  0FF\nM); 
printf("b)  Chromatography  sample  time  is  %.3f  minutes\n", 
usr->runtime) ; 

printf("c)  Retention  time  of  internal  standards\n"); 
for(i  *  1;  i  <*  usr->nstd;  i++) 

printf ("\tstandard  #%d  at  %.3f\n",  i,  (usr->stdrt[i])/convrt) ; 
printf("d)  Areas  of  internal  standards  (used  for  match  calculation)\n"); 
for(i  «  1;  i  <«  usronstd;  i+-*-) 
printf ("\tstandard  #%d  at  %d\n‘‘,  i,  usr->stdarea[i]) ; 
printf ("e)  Minimum  area  required  for  standards  %d\n",  usr-xninarea) ; 
printf ("t)  Change  target  compound  parameters\n" ) ; 
printf ("x)  Return  to  main  menu\n"); 


09-11-87  08:45:46  pmres.c 
Fri  09-11-87  10:04:51 


pmres 


^  #****★*******★★★*★★★★*★■**★•*★★★★*★★★★★**★*■**★*■*-*★*★**★★★★★*★★★★★★**★★★*  **  j 

/*  pmres.c  */ 

/*  send  sample  results  to  printer  */ 

/*  created  6/25/86  B.  Lentz  */ 

/*  V 

^***********************M*************************************************  j 

#include  <stdio.h> 

#  include  “parm.h" 

pmres(pmptr,  position) 

FILE  *pmptr; 
int  position; 

-{ 

extern  struct  value  usr;  /*  run  parameters  */ 
fprintf(pmptr,  "\n\n\n"); 

fprintf(pmptr,  “\n\tSAMPLE  POSITION  #%d\n",  positional); 
swi  tch  (us  r . type [pos i t i on] ) 


case  BLANK: 

— { 

fprintf(pmptr,  "\n\tBlank  "); 
break ; 

case  STAN0RD1 : 

r{ 

fprintf(pmptr,  "\n\tStandard#A  '); 
break; 

L-} 

case  STANDRD2 : 

— { 

fprintf(pmptr,  "\n\tStandard#B  "); 
break; 

L} 

case  STAN0RD3: 

— ( 

fprintf(pmptr,  "\n\tStandard#C  "); 
break ; 

L_} 

case  SAMPLE: 
r-{ 

fprintf (pmptr,  "\n\tWater  Sample  #  %s  ",  usr.nuraber[position]) ; 

break ; 

case  SPIKE: 

r{ 

fprintf  (pmptr,  "\n\tSpike  "); 
break  ; 

L> 

case  SPIKDUP: 


Pg 
of 

1-55 


fprintf  (pmptr,  "\n\t0uplicate  Spike  ”); 
break; 


57 


S?  £ 
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56  default: 

57  r{ 

58  fprintf (pmptr,  '\n\tUnkncwn_Type  ') ; 

59  break; 

60  I—} 

61  — }  /*  end  switch  */ 

62  fprintf  (pmptr,  "\n\tAnalyst  %s",  usr. analyst) ; 

63  fprintf  (pmptr,  u\n\t0ate  %s\t\tTime  %s\n",  usr. date,  usr. time); 

64  fflush(pmptr) ; 


t 


» 


! 


i 


i 

i 


«r.  .w.  -- V\‘  J 
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y*********** *************************************************************/ 


/*  pmseq.c  */ 
/*  print  sample  sequence  to  screen  or  printer  */ 
/*  created  7/24/86  B.  Lentz  */ 

/*  V 


j ******************* *********»******»********•**«*«***********«**•*»*****/ 

♦include  <5tdio.h> 

♦include  "pam.h" 

pmseq(filename) 

char  ‘filename;  /*  device  name  for  redirection  of  output  */ 

— { 

extern  struct  value  usr; 

extern  char  label Q;  /*  date  and  time  information  used  in  name  */ 
int  pos;  /*  sample  position  on  purge  and  trap  */ 

FILE  *fptr; 

if ( (fptr  ■  fopen(filename,  "w"))  *•  0) 

fprintf(stderr,  "\nError  sending  information  to  output  device\n"); 
return; 

L} 

fprintf(fptr,  "\n\n\n\t\tASSIGNED  SAMPLE  NUMBERS'); 
fprintf(fptr,  "\n\tSamples  MUST  be  in  the  order  designed  \n\n"); 
for (pos  -  0;  pos  <  usr.sampcnt;  pos++) 

— ^ 

fprintf(fptr,  "\n\tPosition  %d\t" ,  pos+1); 
switdi(usr.type[pos] ) 

— ( 

case  BLANK: 

r{ 

sprintf (usr. number [pos],  "%d%s.blk",  pos,  label); 

fprintf (fptr ,  "BLANK\n" ) ; 

break; 

L) 

case  STANDRD1: 

— { 

sprintf(usr.nui*)er[pos],  "%d%s.stA",  pos,  label); 

fprintf (fptr,  "STANDARD  A\n") ; 

break; 

*-} 

case  STANDRD2: 

— { 

sprintf(u$r.ruit>er[pos],  "%J%s.stB",  pos,  label); 

fprintf (fptr,  "STANDARD  B\n"); 

break; 

case  STANDR03: 

— { 

sprintf(usr.mmber[pos],  "VJ%s.stC",  pos,  label); 

fprintf(fptr,  "STAN0ARD  C\n") ; 

break; 

case  SAMPLE: 

59 
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56  r{ 

57  fprintf(fptr,  WATER  SAMPLE\t  %s\n“,  usr.number[pos]) ; 

58  break ; 

59  l—} 

60  case  SPIKE; 

61  p{ 

62  sprintf(usr.number[pos] ,  "%d%s.spk'' ,  pos.  label); 

63  fprintf(fptr,  "SPIKE  of  sample  %s\n",  usr.nirtber[pos-l]) ; 

64  break ; 

65  I—} 

66  case  SP1KDUP: 

67  p{ 

68  sprintf(usr.nun*)er[pos3 ,  "%d%s.dsp",  pos,  label); 

69  fprintf(fptr,  "Duplicate  SPIKE  of  sample  %s\n",  usr.number[pos-2]); 

70  break ; 

71  L_} 

72  default: 

73 

74  fprintf(fptr,  "Unasigned  position\n") ; 

75  break; 

76  I—} 

77  L} 

78  L_} 

79  fclose(fptr) ; 

80  I—} 

81 


60 


09-11-87  06:46:14  pmtype.c  Pq 

09-11-87  10:04:51  pnntype  of  1 

1-15 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 


/*  pmtype.c  */ 
/*  print  sample  types  available  to  the  screen  */ 
/*  created  8/13/86  B.  Lentz  */ 

/*  V 


j  **************  *************************************  *************** ******/ 

♦include  <stdio.h> 

pmtypeO 
— { 

printf("\n\t\t\tSAMPLE  TYPE  NUMBERS\n\n\t  0)  BLANK\t\t\t4)  WATER  $AMPLE\n") ; 
printf("\t  1)  STANDARD  #A\t\t5)  SPIKE\n  ) ; 

printf("\t  2)  STANDARD  #B\t\t6)  END  OF  SAMPLES\n\t  3)  STANDARD  #C\n\n"); 


8  5S 
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/*★*★***★**★*******★****+★★*****★*★***★★★**★*********★★*★**********★★*★**/ 


/*  qcexec.c  */ 
/*  calls  QC  data  retrieval  routines  */ 
/*  created  7/18/86  B.  Lentz  V 
/*  V 


j ********************* ****************************** a******************** j 

linclude  <stdio.h> 
main() 

r{ 

printf("\n\nlnsert  QC  Archive  diskette  into  Drive  B\nReady?  (y  or  n):  "); 
whi1e(!reply()) 

printf("QC  Archive  Diskette  in  Drive  B?  (y  or  n):  "); 
printf ("\nUpdating  QC  Archive  Diskette.  Please  wait..."); 
system("copy  a:qc????.dat  b:  >  NUL"); 

printf("\n\nThis  program  will  retrieve  QC  spike  concentrations  \n"); 
printf ("from  the  month  and  year  selected\n"); 
printf ("date  format  month/year  (eg.  7/86)\n"); 
do 

r{ 

qcretrv(); 

printf ("\nDo  you  wish  to  retrieve  more  QC  data?  (y  or  n):  "); 

L-}  while(replyO); 

printf ("\n\nRemove  QC  Archive  diskette  from  Drive  B\nReady?  (y  or  n):  "); 
while(!reply()) 

printf ("QC  Archive  Diskette  removed  from  Orive  B?  (y  or  n):  "); 


09-11-87  08:12:50  calib.c 
Fri  09-11-87  10:04:51  calib 


j *************************** ******* ********** ****** ****** *+**+*++*+*****+ j 


2  /*  calib.c  */ 

3  /*  calib.c  Used  to  calibrate  chromatographic  response  based  on  the  */ 

4  /*  ratio  of  the  response  of  the  target  compound  to  the  standard.  All  */ 

5  /*  three  standards  are  used  in  the  final  concentration  calculations  in  */ 

6  /*  order  to  improve  the  results  (the  three  values  are  averaged).  */ 

7  /*  created  7/24/85  */ 

8  /*  */ 

9  /*  revised  6/13/86  set  calibration  date  */ 

10  /*  B.  Lentz  V 


1 1  /  ***************************  j 

12 

13  finclude  <3tdio.h> 

14  #include  "parm.h" 

15  linclude  "targ.h" 

16  fdefine  CORK  0.90  /*  minimum  required  correlation  */ 

17 

18  calib(uiknwn) 

19  int  unknwn; 
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extern  struct  sample  sampQ 

extern  struct  value  usr; 

int  i,  j,  k; 

float  xmean; 

float  ymean; 

float  xresid; 

float  yresidpWXSTDS]; 

float  ratio  [JUSTUS]  ; 

float  num.- 

float  den; 

float  ycorr; 

float  varl; 

float  r; 

float  intercp; 

float  slope; 

double  sqrt(); 

int  gudcalib; 


j;  /*  pointer  to  target  structure  */ 

/*  set  parameters  */ 

/*  counters  for  loops  */ 

/*  mean  of  x  values  */ 

/*  mean  of  y  values  */ 

/*  x  residuals  V 
/*  y  residuals  */ 

/*  area  ratios  */ 

/*  numerator  */ 

/*  denominator  */ 

/*  sun  of  squares  of  y  residuals  V 
/*  intermediate  value  of  corn  coef  */ 
/*  correlation  coefficient  */ 

/*  calculated  intercept  */ 

/*  calculated  slope  V 
/*  square  root  fuiction  */ 

/*  flag  good  calibration  */ 


ymean  ■  0.0; 
gudcalib  •  TRUE; 

for(k  »  1;  k  <■  usr.nlevels;  k++) 
ymean  +«  samp [unknwn], cone [k] ; 
ymean  •  (float)  ymean/(usr.nlevels) ; 
for(k  •  1;  k  <-  usr.nlevels;  k+-*-) 

— ( 

yresid[k]  •  (float)  samp [unknwn] .cone [k]  -  ymean; 
ycorr  +•  yresidfk]  *  yresid[k]; 

L.} 

for  (j  *  1;  j  <■  usr.nstd;  j++) 


xmean  *  0,0; 

for(k  ■  1;  k  o  usr.nlevels;  k++) 


— { 

ratiofk] 
xmean  +» 


*  (float)  samp(unknwn].area[k]/usr.sarea(k][j); 
ratio[k]  ; 
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1 — } 

xmean  -  xmean/(usr.nlevels); 
mm  *  0.0; 
den  *  0.0: 

for(k  *  1;  k  o  usr.nlevels;  k++) 

r{ 

xresid  *  ratio[k]  -  xmean; 
num  +*  xresid  *  yresidfk]; 
den  +«  xresid  *  xresid; 

slope  -  num/den; 

intercp  *  yroean  -  slope  *  xmean; 
varl  =  sqrt(den/ycorr) ; 
r  *  varl  *  slope; 
if(usr. verbose) 

-{ 

printf("\n  correlation  coeficient  ■  %f\n",  r); 
printf("slope  %f\tintercept  %f\n",  slope,  intercp)' 
-} 

if (r  <  C0RR) 

-{ 

gudcalib  =■  FALSE; 

-} 

else 

-{ 

Scwp[unknwn].$lope[j]  =  slope; 
samp [unknwn] . i ntercp [ j ]  *  intercp; 

-} 


ret urn (gudcalib); 
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I  ******************* **************************************** *************/ 


/*  qcretrv. c  */ 
/*  retrieve  spike  data  from  floppy  disk  files  */ 
/*  created  7/17/86  B.  Lentz  */ 
/*  V 
/*  printed  */ 


j *********** ******** ******iH*******************iHr*************************/ 


# include  <stdio.h> 

♦include  "smonth.h 

qcretrv() 

r{ 

int  flag; 

char  filename[30]; 
int  sday; 
int  eday; 
int  day; 

FILE  *fptr; 

FILE  *pmptr; 
float  cone; 
int  month; 
int  year; 

do 
— ( 

printf("\n\tEnter  month/year  (or  x  to  exit):  "); 
if ((flag  •  readdate(filename,  Smonth,  Syear))  »*  0) 
printf("\nlllegal  date\tRe-enter  month/year:"); 
if (flag  —  -1) 
return; 

— }  while(!valfile(fi1ename)); 
do 
— { 

printf("\nEnter  first  analysis  date  to  view  (day  of  month):  "); 
while((sd*y  ■  readdayO)  ■»  0); 

printf("\nEnter  last  analysis  date  to  view  (day  of  month):  "); 
vdvile((eday  •  readdayO)  •«  0); 

— }  while(sday  >  eday); 

if ( (fptr  »  fopen(filename,  "r"))  ■»  0) 
printf("\n  Error  opening  qc  file\n"); 
else 
— ( 

/*  search  for  starting  date  */ 
do 

if(fscanf(fptr,  "%d  %f",  6day,  Scone)  ««  EOF) 

r{ 

printf("\nHo  QC  data  was  recorded  between  %d  and  %d\n" ,  sday,  eday); 

fclose(fptr); 

return; 

— }  while(day  <  sday); 

if((pmptr  •  fopen("PRN:",  "wb"))  **  0) 

65 


/*  code  to  indicate  type  of  response(valid,  exit)  V 
/*  name  of  qc  file  to  access  V 
/*  starting  day  of  month  to  access  data  */ 

/*  end  day  of  month  to  access  data  */ 

/*  day  of  month  retrieved  from  file  V 
/*  archived  file  */ 

/*  output  to  printer  */ 

/*  concentration  of  spike  sample  V 


Pg  64 
of  105 
1-55 


3 


•A 


3 

•1 

.1 


■j 

•j 

j 

u 

'w 


■V>v. 


-  .  «  /<  .V  .%  .V 


a  a 


i*>u>u>rororv>rorororsjrofv>iVH-tf*—  _ 

N^OOO0^9tU1^wK)MOU}O>NChUlAWNHOU}O9MCriUlAUN^  SS 


qcwrite 


-11-87  06:47:06  qcwrite.c 
-11-87  10:04:51 


/** ********************************** ***************** »*■•  ****************  j 


/*  qcwrite.c 

/*  opens  archive  file  for  writing  spike  concentrations  for  QC  analysis 
/*  qc  filenaaie  format:  qc<J»nthxyear>. dat 
/*  created  6/11/86  B.  Lentz 
/* 


*/ 

V 

V 

V 
*/ 


♦include  «5tdio.h> 

♦include  "para.h" 

qcwrite(conc) 

float  cone;  /*  concentration  of  spiked  sample  */ 

r{ 

extern  struct  value  usr; 

char  fi1ename[20];  /*  name  of  QC  archive  file  */ 

FILE  *fptr,  *fopen(); 

int  day;  /*  day  sanple  was  collected  V 

int  month; 
int  yea*; 

sscanf(usr.date,  "%2d%*c%2<ft*c%2d” ,  Smooth,  Sday,  Syear); 
sprintf (filename,  "a:QC%d%d.dat",  month,  year); 

if((fptr  -  fopen(f11ename.  "a"))  «  0) 

— I 

fprintf(stderr,  "ERROR  opening  control  archive  f11e\n“); 
return; 

L> 

fprintf(fptr,  "%d  %f  ",  day,  cone); 
fclose(fptr); 

L} 
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************************************************************************/ 


/ 

/*  readdate.c 

/♦  wait  for  user  to  enter  in  date  string  and  check  for  valid  entry 
/*  string  should  be  in  format:  month/ year 
/*  created  7/18/86  B.  Lentz 
I* 

^■H,t,*-k*********  **************  ******  *******  ********** 


♦include  <stdio.h> 


readdate(filename,  month,  year) 
char  ‘filename;  /*  name  of  qc  file  */ 

int  ‘month; 
int  ‘year; 

char  string [20] ;  /*  temporary  string  storage  */ 

fgets (string,  15,  stdin); 
if(stmcmp(string,  "x",  1)  *■  0) 
retum(-l);  /*  exit  */ 
sscarf (string,  “%2d%*c%2d" ,  month,  year); 
i f ( ‘month  <  13  &&  *month  >  0) 

sprintf (filename,  b :QC%d«sd .daf  ,  ‘month,  ‘year); 
retum(l)  ; 

return (0) ; 


$ 
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readday. c 


1 


i  y***********************-**********************************#**************/ 


2  /*  readday. c  */ 

3  /*  wait  for  user  input  and  check  for  valid  day  V 

4  /*  created  7/18/86  8.  Lentz  *1 

5  /*  */ 


0  ^**fH*fHH**drtr**dHHHH*******^*****************************flr************l****i»y 

7 

8  linclude  <stdio.h> 

9 

10  readday() 


reply 
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/*  reply. c 

I*  get  yes  or  no  reply  from  user 
/*  returns  true  on  yes  and  false  on  no 
/*  created  7/2/86  8. Lentz 


fdefine  YES  1 
#define  NO  0 
^define  EXIT  -1 


*************************************** i 


V 
*/ 
*/ 

V 

V 

♦******★/ 


reply() 

-{ 

char  string [6] ; 
char  rspns; 
int  value; 


/*  first  letter  of  response  */ 
/*  return  value  */ 


if(!(gets(string,  5))); 
sscanf (string,  "%c",  &rspns); 
if ((rspns  «*  ' y * )  1 1  (rspns  --  * Y* ) ) 
return (YES) ; 

else  if((rspns  *=*  'x')  |j  (rspns  •*  'X')) 
return (EX IT) ; 

else  if ((rspns  -*  'n')  [1  (rspns  -»  'N')) 
return (NO); 
else 

r{ 

printf ("\n\tlllegal  response.  Enter  'y'  or  'n' :  "); 
value  *  'eply(); 
return (value) ; 
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/*  resflpy.c 
/*  store  results  in  floppy  disk  file 
/*  created  6/27/86  B.  Lentz 


7 

7 

7 

7 


/include  <stdio.h> 
/include  "parm.h" 
/include  "targ.h" 


re$flpy(code,  position,  tang,  valuel,  va;ue2) 
int  code;  /*  type  of  results  to  report  7 

int  position;  /*  sample  position  7 

int  targ;  /*  index  of  target  compound  7 

float  valuel,  value2;  /*  concentration  results  7 

r-f 

extern  struct  value  usr; 

extern  struct  sample  sampQ ;  /*  parameters  for  target  cmpds  7 
FILE  *flppy; 

char  fi1ename[20];  /*  floppy  disk  filename  7 
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sprintf (filename,  "b:%s",  usr.nu*er[position]); 
if((f1ppy  -  fopen(filename,  "a"))  0) 

I — { 

fprintf(stderr,  "\nEmor  Opening  %s\n",  filename); 
return; 

M 

switch (code) 

I — { 

case  1:  /*  report  concentration  at  quantitation  limit  7 

I — ( 

fprintf(f1ppy,  "  %s  <  %.3f  %s",  samp(targ].name,  samp[targ].qtlimit,  samp[targ]. units); 
break; 

1 — } 

case  2:  /*  sample  result  7 

I — ( 

fprintf(flppy,  "  %s  ■  %.3f  %s",  sa*p[targ].name,  valuel,  samp[targ]. units); 
break; 

* — } 

case  3:  /*  spike  result  7 


fprintf(f1ppy,  ”  %s  %.3f  %.3f  %s",  samp [targ]. name,  valuel,  value2,  sanp [targ]. units); 

break; 

-> 

case  4: 

I — ( 

fprintf(f1ppy,  “  Recalibration  required"): 

break; 

* — } 

case  5:  /*  standard  calibration  sample  7 

I — ( 

fprintf(f1ppy,  "  Calibration"); 
break; 
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retriev.c 

retrieves  information  from  data  files  on  floppy  disk 


/ 

/ 

/*  created  6/30/86  B.  Lentz 

7* 

/ 

♦include  <stdio.h> 

♦define  MINMATCH  0.7 


retriev(filename) 
char  ‘filename; 

I — ( 

FILE  ‘pmptr; 

FILE  *stptr; 
char  analyst [10]; 
char  date  [10]; 
char  time[10]; 
char  type  [20] ; 
int  npeaks; 
float  rt; 
long  area; 


/*  name  of  data  file  */ 

send  contents  to  printer  */ 
read  values  from  archived  file  */ 
initials  of  analyst  */ 
date  of  analysis  */ 

.  time  of  analysis  */ 

/*  sample  type(blank,  std,  etc)  */ 

/*  number  of  peaks  found  in  sample  */ 
/*  retention  time  in  minutes  */ 
peak  area  */ 


/* 

I* 

/* 

/* 

/* 


/ 


int  i; 

char  string [30];  /*  character  string  V 

char  eguil;  /*  value  is  "<"  or  relationship  */ 

float  valuel,  value2;  /*  floating  point  values  */ 

int  ntargs;  /*  number  of  target  compounds  V 

char  units [10];  /*  units  of  concentration  */ 

;f((pmptr  -  fopen("PRN:",  V*))  -  0) 

'  fprintf(stderr,  "\nError  transmitting  data  to  printer  \n") ; 

return; 

*— } 

fprintf(pmptr,  "\n\n\t%s\n" ,  filename); 
sprintf (string,  "b:%s",  filename); 
i f ( (stptr  «  fopen(string,  "r"))  “  0) 

fprintf(stderr,  "\nError  Opening  %s\n",  filename); 
return; 

L-} 

fscanf(stptr,  "%s  %s  %s  ",  analyst,  date,  time); 

fprintf(pmptr,  "\n\tAnalyst  is  %s\n\t0ate  %s\n\tTime  %s\n",  analyst,  date,  time); 

fscanf (stptr,  "%s  %d",  type,  Snpeaks); 

fprintf(pmptr,  "\tSample  type  %s\n",  type); 

fprintf(prnptr,  "\n  Peak  f\t  retention  time\tarea\n") ; 

for(1  «  0;  1  <  npeaks;  i++) 

i — { 

fscanf (stptr,  "%f  V)",  &rt,  Jarea); 
fprintffpmptr,  "\t%d\t\t%.3At\t%0\n  ,  i‘l.  rt,  area); 

* — ^ 

fscanf (stptr,  "%f  %d",  Svaluel,  Intargs); 
fprintf(pmptr,  "\n\tMatch  factor  is  %.3f\n",  valuel); 
for(i  ■  0;  1  <  ntargs;  i++) 

M  73 
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if(strcmp("B1ank type)  ■«  0) 

rf 

break ; 

M 

else  if(stmcmp("Standard",  type,  8)  **  0) 

p{ 

fscanf (stptr,  "%s",  string); 
if  (stmcnp(string,  "calibration",  6)  •*  0) 


fprintf(pmptr,  "\tCal ibration  run\n"); 


fscanf (stptr,  "%f  %s",  &valuel ,  units); 

fprintf(pmptr,  "\tConcentration  of  %s  is  %.3f  %s\n",  string,  valuel,  units); 

-} 
break ; 


else  if(strcmp("Sample" ,  type)  *■  0) 

-{ 

fscanf(stptr,  "%s  %ls  %f  %s",  string,  Jequil,  Svaluel,  units); 
if(valuel  **  -1.0) 
i  ■  ntargs; 
else 

fprintf(pmptr,  "\tConcentration  of  %s  %c  %.3f  %s\n",  string,  equil,  valuel,  units); 
break ; 

else  if(stmanp("Spik",  type,  4)  -»  0) 


fscanf(stptr,  "%s  %f  %f  %s",  string,  Svaluel,  &value2,  units); 
fprintf (pmptr,  "\tConcentration  of  %s  is  %.3f  %s\n",  string,  valuel,  units); 
fprintf(pmptr,  "\tAmount  of  spike  is.3f  %s\n",  value2,  units); 
break ; 

} 


fprintf  (pmptr,  "\n"); 
wnile(fscanf (stptr,  "%s",  string)  !«  EOF) 
fprintf  (pmptr,  "%s  ",  string); 
fprintf  (pmptr,  "\n\n\n\n")  ; 
fclose(stptr); 
fclose(pmptr) ; 

-} 
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/*******•***********+*****»**********************♦**********»*************  j 

/*  sampler. c  */ 
/*  From  sample  type(std,  blank,  sample),  performs  */ 
/*  necessary  calculations  and  sends  results  to  appropriate  location.  */ 
/*  created  7/24/85  B.  Lentz  V 
/*  V 
/*  revised  5/21/86  added  spike  duplicate  */ 
/*  revised  6/13/86  included  quantitation  limit  */ 
/*  revised  8/21/86  added  dilution  limit  */ 
/*  revised  9/8/86  added  factor  to  indicate  system  response.  Actually  */ 
/*  ratio  of  actual  std  area  to  value  set  in  table  */ 
/*****•**********•*******************************************************/ 
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linclude  <stdio.h> 

♦include  "parm.h" 
linclude  "targ.h" 

sampler(pmptr,  position,  rettm,  area,  calrun) 


FILE  *pmptr;  /*  outp 

int  position;  /*  samp 

int  rettmQ;  /*  sele 

long  areaQ;  /*  sele 

int  *calrun;  /*  cali 

-{ 

extern  struct  value  usr; 
extern  struct  sample  samp Q; 
int  i,  j; 

int  index;  / 

float  ratio;  / 

float  conct[HAXUNKS];  / 

float  sumconc;  / 

float  spikamt;  / 

float  diff;  / 

float  factor;  / 

float  stddev();  / 

float  quant();  / 


/*  output  to  printer  */ 

/*  sample  position  */ 

/*  selected  retention  times  for  std  &  target  peaks  */ 
/*  selected  areas  for  std  &  target  peaks  */ 

/*  calibration  run  */ 


/*  run  parameters  */ 

/*  parameters  for  target  crapds  V 


/*  index  value  to  retrieve  target  info  */ 

/*  ratio  of  target  to  int  std  */ 

/*  concentration  of  unknown  */ 

/*  sum  of  cone  used  for  averaging  */ 

/*  amount  of  target  due  to  spike  */ 

/*  difference  between  expected  and  experimental  V 
/*  indication  of  system  response  */ 

/*  calculates  std  dev  and  returns  mean  V 
/*  calculate  sample  concentration  */ 


swi tch(usr. type [position] ) 

— { 

case  BLANK: 

— { 

/*  check  for  contamination  */ 
for(j  ■  1;  j  <•  usr.nunks;  j++) 
r-{ 

conct[j]  »  quant(j,  area); 

/*  print  out  results  */ 
if (cooct [j]  samp[j].qtlifflit) 

r{ 

/*  contaminated  blank  ♦/ 

fprintf(pmptr,  "\n\t**4  Contaminated  blank  ***\n  ); 
fprintf(stderr,  "\n\t***  Contaminated  blank  ***\n ') ; 
resflpy(6,  position,  j,  conct[j],  0.0); 
exit(0); 

— } 

— } 


I 
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L  break ; 
1 


case  STAND RD1 
case  STANDRD2 
case  STANDRD3 


if(*calrun) 

— { 

fprintf(pmptr,  "\n\tCalibration  run\n"); 
for(i  *  1;  i  <«  usr.nstd;  i++) 

r{ 

usr.sumrtfi]  +=  rettm[i]; 
usr.sumarea[i]  +»  area[i]; 
usr.sarea[usr.type[position]] [i]  *  area[i ] ; 

L-} 

for(i  »  1;  i  <*  usr.nunks;  i++) 

— { 

index  *  usr.nstd  +  i; 

samp[i]  .sumurt  +*  rettm[index] ; 

samp[i].area[usr.type[position]]  *  area[index]; 

if(position  **  3) 
r-{ 

if(usr. invalid) 

— { 

fprintf (pmptr,  "\nRecalibration  is  REQUIRED\n"); 
fprintf(stderr,  "\nRecalibration  is  REQUIRED\n"); 
break; 

>-} 

for(j  *  1;  j  <«  usr.nstd;  j++) 
r-{ 

usr.stdrt[j]  *  usr.sunrt[j]/usr.nlevels; 
usr.stdarea[j]  *  usr.sunarea[j]/usr.nlevels; 

/*  calculate  system  response  factor  */ 
factor*  (float)  area[j]/usr.stdarea[j] ; 
fprintf  (pmptr,  "\tColum  factor  %d  %f\n",  j,  factor); 

M 

for(j  *  1;  j  <*  usr.nunks;  j++) 

r—{ 

sanp[j].rt  *  samp[j].sumurt/(usr.nlevels); 
if (cal ib( j ) ) 

— { 

stocalib(j) ; 

printf("\nSuccessfu1  Calibration  for  %s\n",  samp[jj  .name) ; 
fprintf  (pmptr,  "\n\tCalibration  successful  for  %s\n",  samp(j]  .name) ; 

else 

fprintf  (pmptr,  "\n\t“*  WARNING  ***  Calibration  UNSUCCESSFUl  for  ■  n  .  samp[j]  .name)  ; 
fprintf (stderr,  "\n\t“*  WARNING  ***  Calibration  UNSUCCESSFUL  for  ,s\n  ,  samp[jj .name) ; 

M 

*calrun  *  0; 
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re$flpy(5,  position,  1,  0.0,  0.0); 


for(j  ■  1;  j  <«  usr.nunks;  j++) 

— { 

ccnct[j]  •  quant(j,  area); 
if(conct[j]  <  samp{j]  .qtlimit) 

— ( 

fprintf(pmptr,  "\n\tConcentration  of  %s  is  less  than  %.3f  %s\n",  samp[j] -name,  samp[j].qtlimit,  samp[j]. 
units); 

resflpy(l,  position,  j,  0.0,  0.0); 

L> 

else 

— f 

fprintf(pmptr,  "\n\tConcentration  of  %s  is  %.3f  Vs",  samp  [j]  .name,  conct[j],  samp[j]. units)  ; 
resf1py(2,  position,  j,  conct[j],  0.0); 

L-} 

if ((diff  «  conct[j]  -  samp[j].conc[usr.type[position]))  <  0.0) 
diff  -  -diff; 
if(usr.verbose) 

printf("\naccuracy  diff  %f  resid  %f\n",  diff,  diff-samp[j].accurcy); 
if((diff  -  samp(j] -accurcy)  >0.0) 
r— { 

setcalib(); 

fprintf(pmptr,  "\n***  WARNING  ***  System  out  of  calibration,  recalibration  is  required\n"); 
fpnntf(stderr,  "\n***  WARNING  ***  System  out  of  calibration,  recalibration  is  required\n”); 


case  SAMPLE: 

— ( 

for(j  -  1;  j  <■  usr.nunks;  j++) 

p{ 

conct[j]  •  quant(j,  area); 

$amp[j].lastconc  »  conct[jj; 
if(conct[j]  <  samp[j].qtlimit) 

r{ 

fprintf(pmptr,  "\n\tConcentration  of  %s  is  less  than  %.3f  *s\n' ,  samofil. name,  sampfi) . at! fmit .  sampfil. 
units); 

resflpy(l,  position,  j,  0.0,  0.0); 

else 

r{ 

fprintf(pmptr,  "\n\tConcentration  of  Vs  is  V.3f  Vs",  sampfj]  .name,  conct [j] ,  samp[j] .units) ; 
resflpy(2,  position,  j,  conct[j],  0.01; 
if(conct[j]  >  sampD] -cutoff) 

H 

I  fprintf (pmptr,  "\nWARNlNG  sample  concentration  exceeds  maxiitun  allowable  water  limits\n"); 
if (conct[j]  >  sanp[j]. dilution) 


fprintf(pmptr,  "\nWARNING  sample 
I  fprintf(pmptr,  "For  accurate  rest 


pie  results  are  beyond  the  linear  range  ot  the  instnjment\n"); 
results,  dilute  sample  and  run  again\n); 
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if  (usr.calflag) 

— { 

fprintf(pmptr,  "\n***  WARNING  ***  These  results  are  ONLY  estimates\n") ; 
fprintf(pmptr,  "\nRecali  brat  ion  is  required  for  valid  results’^"); 
resflpy(4,  position,  1,  0.0,  0.0); 

L_} 

break  ; 

<-} 

case  SPIKDUP: 
case  SPIKE: 

r{ 

for(j  *  1;  j  <*  usr.nunks;  j++) 

r{ 

conct[j]  *  quant  (j,  area); 
if((spikamt  »  conctfj]  -  sampjj]. lastconc)  <  0.0) 
spikanrt  *  0.0; 

fprintf(pmptr,  "\n\tConcentration  of  %s  is  %.3f  %s",  samp[j].name,  conct[j],  san^ptjj .units); 
fprintf (pmptr,  "\n\tspiked  amount  of  %s  is  %.3f  %s",  samp[j].name,  spikamt,  samp[j] .units) ; 
resflpy(3,  position,  j,  conct(j],  spikamt); 
qcwrite(spikant) ; 
if(usr.type[position]  »*  SPIKDUP) 

r{ 

if ( (diff  *  conct[j]  -  samp[j] .frstspke)  <  0) 
diff  *  -diff; 

if((diff  -  samp[j].precisn)  >0) 

r{ 

fprintf  (pmptr,  r\n***  WARNING  ***  poor  reproducibility,  recalibrate\n") ; 
setcalib() ; 

M 

L_}  /*  end  SPIKDUP  */ 
i f ( (di ff  «  spikamt  -  samp[j] .spike)  <  0) 
diff  *  -diff; 

if  (di  ff  -  samp[j].accurcy  >0) 

— { 

fprintf  (pmptr,  "\n***  WARNING  ***  poor  sample  recovery,  reca!ibrate\n"); 
fprintf (stderr,  "\n***  WARNING  ***  poor  sample  recovery,  recalibrate\n") ; 
setcalib() ; 

sampfj],  frstspke  =  conct[j]; 

L} 

break; 

default: 

(-{ 

fprintf (stderr,  "XnUlegal  sample  type\n"); 
break ; 

-} 


/*  quantitation  step  */ 
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float  quant(j,  area) 

int  j;  /*  target  anpd  index  */ 

long  areaQ ;  /*  selected  peak  areas  *1 

-{ 

extern  struct  value  usr; 

extern  struct  sample  sanpQ; 

float  wan;  /*  sample  concentration  */ 

int  1; 

int  index;  /*  index  value  to  retrieve  target  info  */ 

float  ratio;  /*  ratio  of  target  to  int  std  */ 

float  conct[MAXUNKS] ;  /*  concentration  of  unknown  */ 

float  suaconc;  /*  sun  of  cone  used  for  averaging  */ 

simconc  ■  0.0; 

for(i  -1;  i  <*  usr.nstd;  i++) 

— { 

index  ■  usr.nstd  +  j; 

ratio  *  (float)  area[index]/area[i]; 

if ( (conct [i]  •  ((samp(j].slope[i])  *  ratio)  +  samp[j].intercp[i])  <  0.0) 
conct[1]  »  0.0; 
surname  +•  conct  [1]; 
if(usr.vertose) 

printf ("unknown  %d  \tconc  Vf\n",  j,  conct [’]) ; 

I—} 

mean  «  stddev (usr.nstd,  sumconc,  conct,  usr.vertose); 
return  (mean); 
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/★A********************************************************************** j 


f*  created  4/1/86 


B.  Lentz 


/*  samptype.c  */ 
/*  sets  sample  lumbers  for  all  sample  types  */ 
/*  created  4/1/86  B.  Lentz  */ 
/*  V 
/*  modified  5/21/86  V 
/************************************************************************/ 


linclude  <5tdio.h> 
linclude  ''parm.b" 


samptype(pos,  type) 

int  *pos;  /*  position  number  */ 

int  type;  /*  sample  type  */ 

-{ 

extern  struct  value  usr; 
char  filename[15]; 


int  code; 
char  string [100]; 
FILE  *fptr; 
int  i ; 


/*  error  code  from  getname 


switch(type) 

r-{ 

case  SAMPLE: 

— { 

usr. type [*pos]  *  type; 

pr i ntf ( " \t\t\t ENTER  sample  IDENTIFICATION:  "); 
while((code  *  getname(filename))  !-  0) 

I  if(code  *■  2) 


if(code 

H 


/*  out  of  space  on  floppy  */ 

/*  delete  files  created  on  current  disk  V 
1;  i  <  *pos;  i++) 


for(i  *  1;  i  <  *pos 
-{ 

if(usr.type[i]  ** 
[-{ 

sprintf (string, 
system(string) ; 
— } 


SAMPLE) 

'del  b:%s  >NUL",  usr.number[i]) ; 


printf ("\n\t\tThere  is  no  space  left  on  the  data  diskette  in  Drive  B\n ' 

printf("\t\tPlease  prepare  a  new  data  diskette\n"); 

system("pause"); 

system("fmtdis")  ; 

system("cls"); 

pmtype() ; 

printf ("\n\tPosition  -»d\tRE-ENTER  sample  IDENTIFICATION:  ,  (4pos)+l); 
*-} 

strcpy(usr.mmber[*pos] ,  filename); 
usr.spikcnt++; 

(*pos)-H-; 

if ((usr.spikcnt  >  9)  &&  (*pos  <  8)) 
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type  -  SPIKE; 
samptype(pos,  type); 


case  SPIKE: 


if(usr.type[(*pos)-l]  *■  SAfff’LE) 

r-{  . 

usr.type[*pos]  «  type; 

printf(''\tPosition  %d\tSPIKE  of  sample  %s\n",  (*pos)+l,  usr.number[(*pos)-l)); 
if (*pos  <  9) 

r{ 

(*pos)++; 

usr.type[*pos]  *  SPIKDUP; 

printf("\tPosition  %d\tOupl icate  SPIKE  of  sample  %s\n",  (*pos)+l,  usr.m*ber[(*pps)-2]); 

(*pos)++; 
usr.spikcnt  -  0; 
break; 

<-} 

else 


fprintf(stderr,  "\tSpike  MUST  follow  sample\n"); 
return; 


case  END: 

— { 

usr.sampcnt  «  *pos; 
*pos  ■  11; 
break; 

default: 


usr. type[*pos]  -  type; 

(*pos)++; 

break; 

} 


return; 

-} 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


/**  ********************************** *******★****************************/ 


/*  select. c  */ 
/*  make  a  numerical  menu  selection  */ 
/*  created  7/10/86  B.  Lentz  V 
/*  V 


j  **★■*★★★■****★**•*★*★★*•**★★*★**★★*★**★* *★*★★★** **★*★★★★**** ★*★*★*★** ★■*****★  j 

linclude  <stdio.h> 
linclude  <ctype.h> 

select/) 

r{ 

char  string [20]; 
char  c; 
long  atoi(); 

while/! (gets (string,  20))); 
sscanf (string,  "%c",  Sc); 
if (isdigit(c)) 
retum(atoi  (string))  ; 
else 

retum(-l) ; 

L-} 
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1 

2 
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4 
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7 

8 

9 
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11 
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21 
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25 

26 

27 

28 

29 

30 

31 

32 

33 
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y ************************************ ***+***********************+**+*++**/ 


/*  seqenc.c  V 
/*  Prompts  user  for  sample  mmbers,  colletion  date/ time  and  */ 
/*  analyst.  It  indicates  which  positions  are  to  be  filled  with  V 
/*  standards,  blanks,  and  samples.  This  information  is  stored  in  a  V 
/*  structure  which  is  passed  to  the  other  routines.  */ 
/*  created  4/17/85  B.  Lentz  */ 
/.  V 
/*  modified  3/31/86  to  allow  more  flexibity  in  loading  seqence  */ 
/*  modified  4/30/86  add  spike  samples  every  10th  sample  */ 
/*  modified  6/10/86  */• 
/*  */ 


♦include  <stdio.h> 

♦include  "parm.h" 

sdqenc(frstnun) 

int  frstrun;  /*  is  this  the  first  run  of  the  day?  */ 

— {• 

extern  struct  value  usr; 

FILE  *  fptr; 

int  pos;  /*  physical  sample  position  */ 

int  type;  /*  sample  type  */ 

char  string [20];  /*  temporary  string  storage  */ 

int  nochg;  /*  nunfcer  of  positions  that  cannot  be  changed  */ 

printf(”\nEnter  analyst  identification  :  "); 
while( I (gets (usr. analyst,  10))) 
fprintf(stderr,  "\nENTER  ANALYST  IDENTIFICATION:  '*); 
printf ("Enter  information  as  requested  below:\n"); 
if(fflush(stdin)  <  0) 
printf ("\nbuffer  not  cleared\n'); 
usr.saae>cnt  ■  10;  /*  number  of  samples  in  run  */ 

pmtype();  /*  print  sample  types  */ 
pos  «  0; 

usr. type [pos]  »  BLANK; 

printf ("\tPosition  %d\tBLANK\n" ,  pos+1); 

pos-w-; 


nochg  »  2; 
if (frstrun) 

— { 

usr.type[pos]  =  STAN0RD1; 

printf ("\tPosition  %d\tSTANDARD#A\n", pos+1) ; 

pos++; 

nochg  ■  3; 

L_) 

v4iile(pos  <  10) 

— { 

printf ("\tPosition  %d\tENTER  sample  TYPE  number:  ",  pos+1); 
type  »  select(); 
while(type  <  0  | !  type  >  6) 
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56  printf("\t\t\tENTER  sample  TYPE  nimber:  ",  pos+1); 

57  type  *  select(); 

58  l—} 

59  samptype(&pos,  type); 

60  1—} 

61  if((fptr  *  fopen(targfile[0] ,  "rw"))  »*  0) 

62  r{ 

63  fprintf(stderr,  "\nERR0R  opening  parameter  fi1e\n"); 

64  return; 

65  l—} 

66  itoa(usr.spikcnt,  string); 

67  wrparm("spikcnt",  string,  fptr); 

68  fclose(fptr) ; 

69 

70  fixseq(nochg,  usr.sampcnt) ; 

71  return; 

72  L} 
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^*rt********************fH»************************************************^ 

/*  setcalib.c  V 

/*  set  switches  to  require  recalibration  on  the  next  run  */ 

/*  created  6/24/86  8.  Lentz  */ 

/*  V 

j *★**********************************************************+******+** **y 


♦include  <stdio.h> 
♦include  "pami.h" 


setcalib() 

-{ 

extern  struct  value  usr; 
FILE  *fptr; 


/*  set  switch  for  calibration  run  */ 
if((fptr  «  fopen(targfile[0],  "rw"))  0) 


fprintf(stderr,  "error  opening  parameter  fi1e\n"); 
return; 

-} 

usr.calflag  »  TRUE;  /*  calibration  required  */ 
wrparm("calflag\  "1",  fptr); 
fclose(fptr) ; 

} 
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Fri  09-11-87  10:04:51  setdate 


/♦♦a*********************************************************************  j 


2  /*  setdate.c  */ 

3  /*  sets  date  and  obtains  operator  id  for  run  */ 

4  /*  created  4/4/86  B.  Lentz  */ 

5  /*  */ 

6  /*  modified  6/6/86  V 

7  /*  modified  7/3/86  V 


g  ^***+********************************************************************^ 

9 

10  linclude  <stdio.h> 

11  #include  "pami.h" 

12 

13  char  label [15];  /*  part  of  std  identification  with  date  &  time  */ 

14 

15  setdateQ 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 


extern  struct  value  usr; 
int  ntime[4] ; 
int  ndate[4] ; 

FILE  *fptr; 

sysdate(ndate,  ntime); 

sprintf (usr. date,  "%2d/%02d/%2d" ,  ndate[0],  ndate[l],  ndate[2]); 
if((fptr  »  fopen(targfile[0] ,  "rw"))  **  0) 

(-{ 

fprintf(stderr,  "error  opening  parameter  file\n"); 
return; 

wrparm("date",  usr. date,  fptr); 
fclose(fptr) ; 

sprintf (usr. time,  "%2d:%02d",  ntime[0],  ntime[l]); 
sprintf (label,  "%02d%02d%02d",  ndate[0],  ndatefl],  ntime[0]); 

-} 
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/M***.******************************************************************/ 

/*  setdio.c  */ 

/*  sends  signal  to  Tekmar  to  control  desorb  start(contact  closure)  V 

/*  open  relay  switch  (set  to  l),or  set  to  0  to  close  relay  switch.  */ 

/*  this  state  will  cause  the  Tekmar  to  wait  for  the  switch  closure  to  */ 

/*  begin  desorb  */ 

/*  created  2/7/86  B.  Lentz  V 

/*  */ 

I  irk  ********  **************************  **********************************  **/ 

♦include  <5tdio.h> 

♦include  "bitset.h"  /*  bit  setting  macros  for  data  acquisition  board  */ 
setdio(value,  port) 

unsigned  char  value;  /*  digital  output  */ 

unsigned  char  port;  /*  set  digital  port  (with  or  w/o  trigger)  */ 

unsigned  char  status;  /*  used  for  error  check  */ 

if  ((STAT  REG  &  0x70)  1-  0) 


fprintf(stderr,  "\nFATAL  ERROR-11  legal  status  register  value\n"); 
fprintf(stderr,  "\nStatus  Register  value  is  %o\n",  STAT  REG) ; 
exit(0); 

} 


COt>_REG(CSTOP); 

status  ■  OATA  OUT; 

while(l  (STAT  REG  8  C0MH  WA1T)); 

COm_REG(CCLEAR); 

while(!  (STAT  REG  &  C0MM  WA1T)); 
C0MM_REG(port) ; 

while(STAT_REG  &  WRITE  WAIT) ; 
OATA_IN(DIOPORT); 

while(STAT_REG  8  URITE_WAIT); 
DATA_lN(value) ; 

while(!(STAT_REG  &  C0ft1_UAIT)); 
status  -  STATJREG; 
if (status  &  0x80) 
ioerr(); 

-} 


'.If. r. 
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j *★***★★******★**★*★★****************★**★***********★**#*★** t* 


stddev.c 

Calculate  std  deviation  of  concentrations  and  throw  out 
outliers.  Acceptance  criteria  -  must  be  within  2  SD. 
created  8/8/85  B.  Lentz 


k ***★*★**★/ 

V 

V 


/************************************************************************/ 
fdefine  SO  2  /*  value  must  be  within  2  std  dev  7 


float  stddev (nvalues,  sun,  cone,  verbose) 


int  nvalues; 
float  sum; 
float  conc[]  ; 
int  verbose; 


/*  number  of  concentration  values  */ 
/*  sum  of  concentration  values  7 
/*  concentration  values  */ 


float  mean; 
int  i,  k; 
float  resid; 
float  sqresid; 
float  stdev; 
float  start; 
float  end; 
double  sqrt(); 


7 


average  concentration 
counters  for  loop  7 
residuals  7 

sum  of  squares  of  residuals  7 
standard  deviation  7 

start  of  window  for  acceptable  cone  values  7 
/*  end  of  window  for  acceptable  cone  values  7 
/*  square  root  function  7 


begin: 

mean  =  sunv'((fIoat)  nvalues); 
sqresid  *  0.0; 

for(i  *  1;  i  <*  nvalues;  i++) 

I — { 

resid  ■  conc[i]  -  mean; 
sqresid  +«  resid  *  resid; 

* — } 

stdev  *  sqrt(sqresid/ (nvalues  -  1)); 
start  *  mean  -  (SD  *  stdev); 
end  *  mean  +  (SD  *  stdev); 
for(i  -  1;  i  <«  nvalues;  i++) 

I — { 

if(conc[i]  <  start  ||  conc[i]  >  end) 
I — ( 

if (verbose) 

printf("outl ier  %f\n",  concfi]); 
sum  «  0; 

for(k  »  1;  k  <  i;  k++) 
sum  +*  conc[k] ; 
for(k  *  i;  k  <  nvalues;  k++) 

-{ 

conc[k]  -  conc[k+l]; 
sum  +*  conc[k] ; 


-} 


if(--nvalues  <  3) 
retum(mean); 
goto  begin; 


68 
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^  ^*********** ********************** A**************************************/ 

2  /'  stdseq. c  V 

3  /*  standard  seqence  is  the  mandatory  sample  seqence  for  the  calibration  */ 


created  3/13/86 


B.  Lentz 


modified  5/20/86 


linclude  <stdio.h> 
linclude  "parm.h" 
linclude  "targ.h" 

stdseq() 

-{ 

extern  struct  value  usr; 
extern  struct  sample  samp[]; 

FILE  *fptr; 

int  pos;  /*  sample  position  on  purge  and  trap  */ 

int  type;  /*  sample  type  */ 

int  i; 

char  string [20]; 

printf("\nEnter  analyst  identification  :  "); 
while(! (gets (usr. analyst,  10))) 
fprintf(stderr,  "\nENTER  ANALYST  IDENTIFICATION  "); 
printf("Enter  information  as  requested  below  for  the  calibration  run:\n"); 

/*  clear  out  input  buffer  */ 
i f (ff lush(stdin)  <  0) 
printf("\nbuffer  not  cleared\n"); 
usr.sampcnt  -  10;  /*  number  of  samples  in  run  */ 

pm  type  ()  ; 

usr.type[0]  *  BLANK; 

printf("\tPosition  l\tBLANK\n") ; 

usr.type[l]  -  STANDRD1; 

printf ("\tPosition  2\tSTANDARD#A\n" ) ; 

usr. type [2]  ■  STANDRD2; 

printf ("\tPosition  3\tSTANDARD#B\n" ) ; 

usr. type [3]  «  STANDRD3; 

printf ("\tPosition  4\tSTAN0ARD#C\n") ; 

for(pos  *  4;  pos  <  10;) 

— { 

printf ("\tPosition  %d\tENTER  sample  TYPE  number:  ",  pos+1); 
type  *  selectO; 
while(type  <011  type  >  6) 

— ( 

printf ("\t\t\tENTER  sample  TYPE  nunber:  ,  posM); 
type  *  select'); 

L> 

samptype(&pos,  type); 

i  L-} 

if ((fptr  *  fopen(targfile[0] ,  “rw"))  -*  0) 
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56-73 

56  fprintf(stderr,  "\nERR0R  opening  parameter  file\n"); 

57  exit(O); 

58  L} 

59  itoa(usr.spikcnt,  string); 

60  wrparm("spikcnt",  string,  fptr); 

61  fclose(fptr); 

62 

63  fixseq(5,  usr.sampcnt); 

64  for(i  ■  1;  i  <«  usr.nstd;  i++) 

65 

66 

67 

68 

69  for(i  «  1;  i  <*  usr.nunks;  i++) 

70  samp[i].sumurt  *  0; 

71  return; 

72  L} 

73 


— { 

usr.sumrt[i]  »  0; 
usr.sumarea[i]  *  0; 
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^************* ********************************************** *********** **  j 


/*  stocalib.c  */ 
/*  change  values  in  "permanent"  parameter  file  for  slope  and  intercept  */ 
/*  and  change  calibration  flag  to  indicate  that  calibration  was  */ 
/*  successful  */ 
/*  created  8/6/86  B.  Lentz  */ 

/*  V 


y ******************* ***************** ************************************y 


♦include  <5tdio.h> 

♦include  "parm.h" 

♦include  "targ.h" 

stocalib(irknwn) 
int  unknwn; 

{ 

extern  struct  sample  sanp[];  /*  pointer  to  target  structure  */ 
extern  struct  value  usr;  /*  set  parameters  */ 
int  j  ; 

FILE  *fptr; 

char  string[50];  /*  string  representation  of  number  V 

char  pamrame[2G] ;  /*  parameter  name  */ 

if((fptr  «  fopen(targfi le [unknwn] ,  "rw"))  »•  0) 

— ( 

fprintf(stderr,  "\nERR0R  opening  parameter  file  %s\n",  targfile[irknwn]); 
return; 

for  (j  *  1;  j  <•  usr.nstd;  j++) 

r{ 

ftoa(samp[unknwn].slopetj] ,  string,  4,  4); 
sprintf(pamname)  "s1ope%d",  j); 
wrpamKparmname,  string,  fptr); 
ftoa(sanp[unknwn].intercp[j],  string,  4,  4); 
sprintf(parmname,  "intercp%d",  j); 
wrpanm(parmname,  string,  fptr); 

fclose(fptr) ; 

if  ((fptr  »  fopen(targfile[0],  "rw"))  *»  0) 

fprintf (stderr,  "\nERR0R  opening  parameter  file  %s\n",  targfile[0]); 
return; 

M 

usr.calflag  *  FALSE; 
wrparm("calf1ag",  "0",  fptr); 
strcpy(usr.caldate,  usr.date); 
strcpyfstring,  usr.caldate); 
wrparmC'caldate",  string,  fptr); 
fdose(fptr)  ; 
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/M1H»»!»*IH>*********  ************»*****************•*******************•**/ 

/*  stoflpy.c  */ 

/*  store  sample  info,  retention  time  end  areas  on  floppy  disk  "b" ,  V 
/*  created  6/13/86  B.  Lentz  V 

/*  V 

j  ******************************************************  *******************  j 


finclude  <stdio.h> 
#  include  "parm.h" 


stoflpy(position,  npeaks, 
int  position; 
int  npeaks; 
int  rettm[MAXSTDS] ; 
float  area[MAXSTDS] ; 
float  hit; 

— ( 

extern  struct  value  usr; 
FILE  *stptr; 
char  filename[20]; 
char  strtyp[20] ; 
float  rt; 
int  i; 


rettm,  area,  hit) 

/*  saeple  position  V 

/*  number  of  peaks  found  in  sample  */ 

/*  sample  retention  tines  */ 

/*  saeple  areas  */ 

/*  natch  factor  */ 

/*  run  parameters  */ 

/*  to  send  data  for  archiving  on  floppy  */ 
/*  filename  on  floppy  */ 

/*  string  containing  sample  type  V 
/*  retention  time  in  minutes  */ 


swi tch (usr . type [pos i ti on] ) 
[1*.  BLANK: 


— { 

strcpy(strtyp,  "Blank  "); 
break; 

case  STANDRD1: 

— ( 

strcpy(strtyp,  "Standard# A  "); 
break; 

case  STANDRD2: 

— { 

strcpy(strtyp,  "Standard#B  "); 
break; 

L_} 

case  STAN0RD3 : 

— { 

strcpy(strtyp,  "Standardise  "); 

break; 

M 

case  SAMPLE: 


— ( 

strcpy(strtyp,  "Sample  "); 
break; 

case  SPIKE: 

— ( 

strcpy(strtvp,  "Spike 
break; 
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* — } 

case  SPIKDUP: 

I” strcpy(strtyp,  " Sp i ke_0up  1  i cat e  "); 


r{ 

strcpy 

break; 

*-} 

default: 


r~{ 

strcpy 

break; 

—1 


strcpy(strtyp,  "lkiknown_Type  "); 


I — }  /*  end  switch  */ 

sprintf(filename,  "b:Vs",  usr.number[position3); 
if((stptr  *  fopen (filename,  "w"))  0) 

— { 

fprintf(stderr,  "\nError  Opening  %s\n",  filename); 
return; 

fprintffstptr,  "Vs  Vs  Vs  ",  usr. analyst,  usr.date,  usr.time); 
fprintf(stptr,  "Vs  ",  strtyp); 

fprintf(stptr,  "Vd  ",  npeaks); 
for(i  -  1;  i  <*  npeaks;  i++) 


rt  -  rettm[ij/(usr.freq  *  60.0); 
fprintf(stptr,  "V.2f  V.Of  ",  rt,  area[i]); 


fprintf(stptr,  "V.3f  Vd  ",  hit,  usr.nunks); 
fclose(stptr) ; 

-} 
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f*  sysdate.c 

/*  get  current  system  date  and  time 
/*  created  8/7/85  8.  Lentz 

/* 


$ysdate(date,  time) 
int  dated; 
int  timed; 

struct  regval  {  int  ax.bx.cx.dx.si.di.ds.es;  }  srv; 
int  ndat[4]; 

int  year;  /*  4  digit  listing  of  year  V 

char  string [6];  /*  string  representation  of  year  V 

srv. ax  *  0x2 aOO; 

sysint21(8srv,&srv); 

year  «  srv.cx;  /*  year  */ 

sprintf (string,  "*d",  year); 

sscanf(string,  "%*c%*c%2d" ,  i (date [2])); 

ndat[l]  -  srv.dx; 

srv. ax  «  0x2c00; 
sys int21 (&srv , Ssrv) ; 
ndat[2]  -  srv.cx; 

date[0]  •  ndat[l]  »  8;  /*  mnrtfi  */ 

dateCl]  *  ndat(l]  &  Oxff;  /*  day  */ 

time[0]  •  ndat[2]  »  8;  /*  hour  */ 

time[l]  •  ndat[2]  &  Oxff;  /*  minutes  */ 


*/ 

V 

*/ 

V 
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♦include  <3tdio.h> 

♦include  "pann.h" 

temp(position,  smprt,  areas) 
int  position; 

int  smprt [NPEAKS] ;  /*  retention  times  of  sample  peaks  */ 

float  areas [NPEAKS] ;  /*  areas  of  san^le  peaks  */ 

-{ 

extern  struct  value  usr; 

FILE  *fptr; 

int  npeaks;  /*  nimber  of  peaks  found  in  sample  */ 
int  j; 

float  varl;  /*  retention  time  in  minutes  */ 

— static  char  *file[]  •  { 

"blank.dat", 

"standrdl.dat ", 

"standrd2.dat” , 

"standrd3.dat", 

"sample.dat", 

"samplel.dat", 

"sample2.dat" , 

"sample3.dat", 

"sample4.dat", 

"sample5.dat" 

I—}; 

if((fptr  *  fopen(file[position] ,  "r"))  *■  0) 
r-{ 

fprintf(stderr,  "\nError  opening  %s  datafile\n" ,  file[position]); 
exit(0) ; 

— } 

fscanf(fptr,  "%d  ",  ftnpeaks); 
for(j  *  1;  j  <*  npeaks;  j++) 

fscanf(fptr,  "%f  ",  &varl); 
smprt[j]  ■  varl  *  60.0  *  usr.freq; 
fscanf(fptr,  "%f  ",  Sareas [j]); 

fclose(fptr); 
retum(npeaks) ; 

-}  /*  end  temp  V 


,  t 


main 
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#inc!ude  <stdio.h> 
main() 

r{ 

char  buffer[6]; 

printf("\n\nlnsert  a  new  diskette  into  Orive  B\nPress  Enter  key  when  ready..."); 

fgetc(stdin) ; 

fflush(stdin); 

printf("\nPlease  Wait..."); 

pri  ntf  ("Enter"); 

f gets (buffer,  5,  stdin); 

fflush(stdin); 

pri  ntf  ("Enter"); 

fgetc(stdin); 

L-} 
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j ************************************************************************ j 


User  menu 


this  is  the  main  menu  routine  to  control  sample  automation  on  the  gc  */ 


/*  Created  3/28/85  B.  Lentz  */ 
/*  */ 
/*  last  revision  4/21/86  */ 

/ *HkM***Ht  **************************************************************** j 


8.  Lentz 


♦include  <stdio.h> 
♦define  TRUE  (0--0) 
♦define  FALSE  (TRUE 


main(angc,  argv) 
int  argc; 
char  *argv[]; 

-{ 

int  instmo; 


/*  menu  selection  */ 


trfiile(0**0) 

I  system("c1s"); 


/*  infinite  loop  *1 


system("c1s") ;  /*  clear  screen  */ 

printf("\n\n\n\n\n\n\t\tSelect  one  of  the  following  options:\n"); 

printf ("\n\n\t\t  1)  Analyze  for  TCE  \n") ; 

printf("\t\t  2)  Calibration  run  for  TCE\n"); 

printf ("\t\t  3)  Prepare  new  Data  diskette(for  drive  B)\n“); 

printf("\t\t  4)  Prepare  new  VGA  Data  System  diskette(for  drive  A)\n"); 

printf("\t\t  5)  Retrieve  QC  data\n"); 

printf ("\t\t  6)  Retrieve  archived  data\n"); 

printf ("\t\t  7)  Change  operating  parameters \n" ); 

printf ("\t\t  9)  Exit  Menu\n"); 

printf ("\n\n\t  ENTER  OPTION  NUMBER:\t") ; 

instmo  «  selectQ; 

vrfii  le( instmo  <  1  |[  instmo  >  9) 


fprintf(stderr,  "Invalid  Input.  Enter  NUMBER  of  selected  option:\t"); 
instmo  •  selectQ ; 


system("cls");  /*  clear  screen  */ 

switch  (instmo) 

r{ 

case  1:  /*  sample  run  */ 

r{ 

printf ("Analyze  for  TCE\n" ) ; 
system("analyz  0"); 
break; 

— }  /*  end  case  2  V 

case  2:  /*  calibration  run  */ 

r{ 

printf ("Calibration  run  for  TCE\n" ) ; 
system("analyz  In¬ 
break; 

case  3:  /*  format  data  diskette  */ 

|  printf ("Prepare  data  diskette  for  u$e\n"); 
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76 

77 

78 

79 

80 
81 
82 
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system("fmtdis") ; 
break; 

} 

case  4:  /*  Format  system  floppy  disk  */ 

I _ { 

printf("\nPrepare  new  VOA  Data  System  diskette\n"); 
system("copsys") ; 
break; 

I — } 

case  5:  /*  QC  V 

{ 

printf ("Retrieve  QC  0ata\n"); 
system(  "qcexec") ; 
break; 

1 — } 

case  6:  /*  retrieve  data  from  floppy  */ 

I — { 

printf ("Retrieve  archived  data  from  diskette\n"); 
system(" xretrv") ; 
break; 

I — } 

case  7:  /*  Modify  operating  parameters  V 

r_{ 

printf  ("Change  operating  parameters^" )  ; 
system("change") ; 
break; 

I — } 

case  9:  /*  exit  to  DOS  */ 

I — { 

printf ("\nRetuming  to  operating  system\n"); 
printf("Type  \"menu\"  to  return  to  data  analysis  system\n"); 
exit(0); 
break; 

-} 

default: 
r~{ 

printf("\nNot  a  valid  selection\n"); 

break; 

} 

«-} 

!_} 

-} 
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j *************************★******★**★♦★***★******************************/ 


/*  valfile.c  */ 
/*  check  for  presence  of  file  */ 
/*  created  7/18/86  B.  Lentz  */ 

/*  ^  ^  *****,***.♦*!/ 


♦include  <stdio.h> 

valfile(filename) 
char  ‘filename; 

[-{ 

FILE  *fptr; 

if((fptr  *  fopen(filename,  "r"))  «■  0) 

— { 

printf( "\nFile  %s  does  not  exist\n",  filename); 
return (0); 

L} 

fclose(fptr) ; 
return(l) ; 


'in?  7* .» VL»  77  ^ 
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/*  wrparw.c 

/*  Write  new  value  of  parameter  to  "panu.dat" 
/*  created  8/16/85  B.  lent! 

/* 

/*  wodlfled  6/9/86 


♦include  <std1o.h> 
♦define  BASE  0 


wrparm(nane,  newval,  fileptr) 
char  *na»e;  /*  name  of  parameter  to  change  */ 
char  ‘newval;  /*  new  value  of  par awe  ter  V 
FILE  ‘fileptr;  /*  pointer  to  file  to  change  */ 


char  ‘newval; 
FILE  ‘fileptr; 

-{ 

char  string(20j; 
char  oldval  [20] ; 
long  pos; 

rewind(flleptr); 

do 


par— ter  nawe  frow  ■pem.det"  •/ 


/*  read  par— ter  n— 
/*  *3  par—  */ 

/*  position  In  flit  */ 


if(fscarf (fileptr,  V.  string)  «  0) 

— { 

printf("\p%t  mt  tamd  in  fHaV>*.  String); 
return; 


— }  while(stra*p(string,  n— )  !•  0); 
pos  *  ftell (fileptr); 
fscanf (fileptr.  "%s  %s",  string,  oldval); 
if(fseek(fileptr,  pos.  BASE)  "■  -I) 


prlntf (“Address  error  -  laptop er  swtk\n"); 
return; 


f prlntf (fileptr.  "  \t%-15$  %-15s“,  newval.  oldval); 
if(fflush(fileptr)  ■■  -1) 
fprintf(stderr,  "ERROR  writing  to  f11e\n“); 
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/*  xretrv.c 

/*  executive  module  for  retrieving  data  from  floppy  disk 
/*  created  7/1/86  8.  Lentz 

/* 

/* 

main() 


[r 


filelist(); 

} 


V 

V 

V 

*/ 


of 
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smooth 


I 

f 

f 

e 

* 

p* 

pi 

► 

f 


56 

ctr  *  6; 

57 

break ; 

58 

M 

59 

case  13: 

60 

r{ 

61 

none  *  143; 

62 

a[0]  *  25; 

63 

a[l]  •  24; 

64 

a[2]  *  21; 

65 

a [3]  *  16; 

66 

a  [4]  ■  9; 

67 

a[5]  *  0; 

68 

a  [6]  •  -11; 

69 

ctr  *  7; 

70 

break; 

71 

*— } 

72 

case  15: 

73 

r-{ 

74 

norm  *  1105; 

75 

a[0]  -  167; 

76 

a [1]  *  162; 

77 

a [2]  •  147; 

78 

a [3]  -  122; 

79 

a [4]  *  87; 

80 

a [5]  •  42; 

81 

a [6]  *  -13; 

82 

a [7]  *  -78; 

83 

ctr  *  8; 

84 

break; 

85 

l— } 

86 

case  17: 

87 

r{ 

88 

norm  *  323; 

89 

a[0]  ■  43; 

90 

a[l]  *  42; 

91 

a [2]  ■  39; 

92 

a [3]  •  34; 

93 

a [4]  •  27; 

94 

a [5]  •  18; 

95 

a [6]  *  7; 

96 

a [7]  *  -6; 

97 

a [8]  *  -21; 

98 

ctr  *  9: 

99 

break ; 

100 

L-} 

101 

case  19: 

102 

r{ 

103 

norm  *  2261 

104 

a[0]  *  269; 

105 

a[l]  *  264, 

106 

a  [2]  •  249; 

107 

a [3]  *  224; 

108 

a [4]  •  189; 

109 

a[5]  *  144; 

110 

a (6]  *  89; 

> 
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a [7]  ■  24; 
a  [8]  -  -51; 
a [9]  -  -136; 
ctr  ■  10; 
break; 

case  21: 


norm  «  3059; 
a[0]  «  329; 
a[l]  »  324; 
a [2]  ■  309; 
a [3]  ■  284; 
a [4]  ■  249; 
a [5]  ■  204; 
a [6]  •  149; 
a [7]  -  84; 


a [9]  •  -76; 
a[10]  •  -171; 
ctr  *  11: 


case  23: 

— ( 

non*  -  805 
a[0]  -  79; 
a 

a  [21  ■  75; 


a  [4]  •  63; 


a [6]  •  43; 
a  [7]  ■  30; 


a [10]  »  -21; 
a[ll]  -  -42; 
ctr  *  12; 


norm 

«  5175 

a[0j 

-  467; 

a[l] 

•  462; 

a[2] 

-  447; 

a[3] 

*  422; 

a  [4] 

*  387; 

a  [5] 

•  342; 

a  [6] 

«  287; 

a  [7] 

•  222; 

a  [8] 

•  147; 

« n*  r\*  i\*  -ji  7*r+Lni'jyrcv^r*r*.v‘.  'AW/JVJV^v.VTRWlTOVW^JV^VATOV.Tr.v'^AV"/.v.v.'v,Jv.T,.-v,.v.v.v.v.'A^ 
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a[ll]  ■  -138; 
a [12]  •  -253; 
ctr  *  13; 
break; 

case  9: 
default: 

I — { 

norm  «  231; 
a[0]  •  59; 
a[l]  -  54; 
a [2]  -  39; 
a [3]  -  14; 
a [4]  -  -21; 
ctr  *  5; 
factor  *  9; 
break; 

M 
1 — } 

for(i  «  ctr;  i  <  13;  i++) 
a[i]  ■  0; 

for(i  ■  2;  i  <■  factor;  i++) 
np[i]  «  data[i-2]; 
spts  ■  npts  -  ctr  -  1; 
for(i  •  ctr;  i  <  spts;  i++) 

r~ { 

j  »  i  +  ctr  -  1; 
for(k  *  1;  k  <  factor;  k-H-) 
np[k]  •  np[k+l]; 
np  [factor]  -  data[j]; 

nsum  •  a[0]  *  np[ctr] ; 
for(x  ■  1;  x  <  ctr;  x++) 
nsun  +•  a[x]  *  (np[ctr+x]  +  np[ctr-x]); 
data[i]  *  nsum/norm; 
l_} 
return; 

>-} 
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